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ABSTRACT 


This thesis involved a numerical experiment to compare a _ deterministic 
Generalized Lanchester Equation model, referred to as the M/W model, and a 
stochastic computer simulation model, referred to as the C/S model. A discussion of 
the historical background of Lanchester’s equations precedes the presentation of the 
two models and the experimental design. The results are presented graphically and 
show that the M/W force level trajectory is a good approximation for the C/S force 
level trajectory. It was also shown that the two model's trajectories behaved similarly. 
These results indicate that deterministic attrition models may often be good 
approximations for the mean of stochastic attrition models. Command and control 
applications of a model like the M/W model, are presented and a list of suggested 


follow-on research is provided to stimulate further work in this area. 


THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may 
not have been exercised for all cases of interest. While every effort has been made, 
within the time available, to ensure that the programs are free of computational and 
logic errors, they cannot be considered validated. Any application of these programs 


without additional verification is at the risk of the user. 
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I. INTRODUCTION 


A. PURPOSE AND OVERVIEW 

The main purpose of this thesis is to conduct an experiment, using two different 
types of combat attrition models, to determine whether their results are similar. There 
is a belief that deterministic Lanchester Equation type models closely parallel the mean 
results of stochastic computer simulations. This research will attempt to determine 
whether the results of a deterministic “Generalized Lanchester Equation” type model 
and a stochastic computer simulation model are similar. This work is the first such 
research of this type published to the author’s knowledge. 

Chapter One provides the reader with an introduction to modeling and an 
overview of the history of Lanchester Equation type models. The purpose of Chapter 
Two 1s to introduce the two models which are utilized in this experrment and Chapter 
Three is a presentation of the experimental design. The results of the experiment are 
provided in Chapter Four and a discussion of the conclusions and recommended 


follow-up research are given in Chapter Five. 


B. MODELING 
Throughout history, military analysts have been developing and utilizing various 
types of models for detailed analysis to assist in decision making. A model is defined 
as a “simplified representation of the entity it imitates or simulates” [Ref. 1: p. 1]. A 
second definition of a model used by the U.S. Army Models Review Committee, “an 
abstract representation of reality which is used for the purpose of prediction and to 
develop understanding about the real-world processes” [Ref. 2: p. 5], implies that 
models are designed to be close representations of a real-world entity or process. A 
final definition for a military model is, “an abstraction of reality, the elements of which 
are chosen for (a) an investigative purpose or (b) a resource management purpose; in 
other words, an abstraction to assist in making decisions.” [Ref. 1: p. 3]. From these 
definitions, it can be seen how important the development and use of models is to the 
military analyst. 
In general, there are three types of models: 
1) Iconic model 
2) Analogue model 
3) Symbolic model 


An iconic model is a miniature version of the entity, such as an airplane or a tank 
model. An analogue model is an artificial representation of reality. An example of an 
analogue model is a map which represents the three dimensional real-world on a two 
dimensional, small scale sheet of paper. A symbolic model is one in which words or 
numerical descriptions are used to represent an entity or process. An example of a 
symbolic model is a mathematical equation which represents a process such as 
attrition. 

The types of models used by mulitary analysts cover the entire spectrum from 
military field exercises (Which can be thought of as iconic models) to analytical models. 
Figure 1.1 is a combination of information found in Tavlor’s book Force-on-Force 
Aurition Modelling {Ref. 2: p. 7], and Hughes’ book Military Modeling {Ref. 1: p. 10}, 
which shows the various tvpes of models in use and their characteristics of operational 


realism, degree of abstraction, convenience and accessibility. 


cel cll enn le nn nella eeeneeeaetereerreeremernmnmamenliith 
MILITARY MILITARY MAP WAR COMPUTER ANALYTICAL 

FIELD FIELD EXERCISES GAMES SIMULATIONS MODELS 
EXERCISES EXPERIMENTS 


HUMAN DECISION IMPACT AND OPERATIONAL REALISM 


DEGREE OF ABSTRACTION 


OUTCOME REPRODUCIBILITY 


CONVENIENCE AND FLEXIBILITY 





Figure 1.1 Current Model Types. 
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The three most convenient and accessible types of models shown in Figure 1.1 
are War games, computer simulations and analytical models. These three types of 
models are also the most commonly used to model combat attrition. A brief 
description of each of them will be given to better understand and help distinguish 
between them. 

1. War Games 

War games are conducted using real-world combat scenarios and people 
playing various positions on a headquarters battle staff. A war game is conducted to 
allow a commander and his staff to exercise the staff decision making processes in 
‘realistic’ combat scenarios. This allows the commander to exercise and develop his 
staff without incurring the true cost of combat, the loss of human life. War games can 
be run utilizing the sand table method or on computers. The key to war gaming is that 
there is always a person making the decisions that are required in the game. 

2. Computer Simulations 

Computer simulations are also used to model the combat attrition process. 
The primary difference between computer simulations and war games is the method 
used to represent the decision process. War games use people playing staff roles to 
make the decisions where simulations use algorithms to represent the decision process. 
A combat simulation not only represents the combat process, it acts it out from start 
to finish. A computer simulation begins with a set of input parameters and runs 
continuously until completion of the battle and then provides the results of the battle 
as output. Computer simulations which utilize pseudo-random number generators to 
determine the results of random events, such as the outcome of one soldier firing at 
another, are called Monte Carlo Simulations and are stochastic models. 

3. Analytical Models 

Analytical models are the third type of model used in modeling combat 
attrition. Analytical models are symbolic models which use mathematical symbols and 
equations to represent the combat attrition process. Analytical models can be 
developed as stochastic models or deterministic models. Stochastic models, as 
introduced in the previous section, utilize probablity distributions to determine certain 
variables. Therefore the output from two consecutive runs of a stochastic model will 
more than likely be different. The deterministic model however, will produce the same 


set of output values for a given set of input parameters. 


14 


My research 1s focusing on the comparison of a stochastic Monte Carlo 
Simulation and a deternunistic analytical model. The deterministic model | am using 
for comparison is a Generalized Lanchester Equation model which was developed from 
the original Lanchester Equation Theory. To provide a basis for the readers 
understanding of this experiment, I will give a brief introduction to the history of the 


Lanchester Equation type models. 


C. LANCHESTER EQUATION THEORY 
1. History Of Lanchester Equation Type Modeling 

F. W. Lanchester was a British aeronautical engineer, who in 1914 developed 
two sets of simple differential equations to model the combat attrition process between 
two opposing homogeneous forces (1.e., fighter aircraft vs.. fighter aircraft). His 
intentions were to provide insight into the dynamics of combat under ‘modern 
conditions’ of warfare and to justify the principle of concentration of forces. The term 
modern conditions of warfare will be addressed again later. 

Lanchester’s original work was designed to model a force-on-force attrition 
process involving two homogeneous forces. He developed two combat attrition 
models, each having it's own unique assumptions, for this purpose. These original 
models were designed to model an ‘aimed fire’ combat scenario and an ‘area fire’ 
combat scenario. These two models are often referred to as the classical Lanchester 
Equation Theory. 

2. Modeling Aimed Fire 

The phrase aimed fire, as Lanchester used it in his original work, refers to 
combat between two forces, X and Y, where each combatant from X force acquires a 
Y force target (1.e., locates and takes aim) and fires. An example of this type combat is 
an infantry battle for control of an area such as hill 224. The following assumptions 
apply for this model: 

ieee Lhe combat being modeled involves two homogeneous forces (i.e., Infantry vs. 
Infantry). 

2) The entire X force and Y force are within weapons range of one another. 

3) The effects of weapons rounds are independent. 

4) Each of the forces is well enough aware of the location and condition of all 
enemy forces so that they will engage only live enemy units. Also the rate at 
which they kill enemy targets is constant. 


5) Fire is uniformly distributed over surviving enemy targets. 


[> 


Lanchester’s original model for aimed fire is given by Equation Set 1.1: 


i 


dx/dt = -ay(t) with x(0)= X,, 


(eqnule) 


dy/dt = -bx(t) with y(0)= Yo5; 

where x(t) is defined as the number of combatants of X force at time t, y(t) 1s defined 
as the number of combatants of Y force at time t, and a and b are constants which are 
called attrition-rate coefficients. The attrition-rate coefficient ‘a’ represents the 
effectiveness of a Y combatant killing a X combatant per unit of time. Similarly the 
attrition-rate coefficient ‘b’ represents the effectiveness of a X combatant killing a Y 
combatant per unit of time. This set of equations assumes that the rate at which a 
force X can cause casualties to a Y force 1s proportional to the number of combatants 
in the X force. Similarly, the rate at which a combatant from Y force can cause 
casualties to the X force is assumed to be proportional to the number of combatants in 
the Y force. When Equation Set I.l_1s integrated, the resulting state equation is 


obtained and has been labeled the Lanchester Square Law: 
b(X9-x(t)7) = alyo7-V(t)?). (eqn 1.2) 


3. Modeling Area Fire 
The phrase area fire as Lanchester used it in his original work, refers to a 
combat scenario involving two homogeneous forces, say X and Y, which are uncertain 
of exact enemy locations and engage one another by firing in the general area where 
the enemy force is located. This scenario has the following assumption set: 

1) The combat being modeled involves two homogeneous forces (1.e., Infantry vs. 
Infantry) 

2) The entire X force and Y force are within weapons range of one another. 

3) The effects of weapons rounds are independent. 

4) Each force is aware only of the general location of the enemy force and 
therefore engages the enemy by firing into that general ‘area’ without the 
benefit of knowing their effectiveness. 

5) Fire from all surviving combatants is uniformly distributed over the area 
which the enemy occupies. 


6) Each force has the same vulnerable area to enemy fire. 


16 


Lanchesters original model for area fire is defined by the Equation Set 1.3: 


dx/dt = -axy, 


(eqn 1.3) 


dy/dt = -byx. 
This model assumes that the rate of attrition is proportional to both force levels x(t) 
and y(t). This means that in the area fire model, the rate at which a X force can cause 
casualties to a Y force is not only dependant upon the number of combatants in its 
own force, but also on the number of combatants alive in the Y force. The rate at 
which the Y force can cause casualties on the X force is similarly dependent on both 
their own force level and the force level of the opposing force X. 

The well Known ‘Linear Law’ is obtained by integration of Equation Set 1.3 


and is shown here in Equation 1.4: 


B(X_-X(t)) = aly,-y(t)), (eqn 1.4) 


where X, and y, are the initial force levels at time t=O, x(t) and y(t) are the force 
levels at time t, and a and b are the attrition-rate coefficients for each force. 
4. Attrition Rate Coefficients 

The Lanchester attrition-rate coefhicients were introduced in the models above 
and were defined as constants. In order to use the Lanchester models there must exist 
a set of acceptable attrition-rate coefhcients for each possible type of engagement. 
Lanchester’s original model utilizing homogeneous forces, only requires the availability 
of a limited number of coefficients. However, in Lanchester Equation type models 
which allow for the combat between heterogeneous forces, the number of attrition-rate 
coefficients grows very rapidly when one considers all the possible combinations of 
combatant/weapons vs. combatant/weapons that can occur. As Lanchester postulated 
his original work, he used constants for attrition-rate coefficients, which implies that 
the kill capability of a force does not change over time. With constant attrition-rate 
coefficients these models are easy to solve. However, there is no reason why we should 
believe that such factors as range, weather, visibility, and training do not affect the kill 
capability of a combatant. Taylor provides two methods for determining numerical 
values for Lanchester attrition-rate coefficients which are used in the United States. 
[Ref: 2: p. 45] 
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1) A Statistical estimate based on ‘combat data generated by a detailed Monte 

Carlo combat simulation. 

2) An analytical submodel of the attrition process for the particular combination 
of firer and targets. 

In the first method, the simulation output is used to fit one or more free 
parameters in the analytical model in an attempt to have it to provide similar results. 
S. Bonder has labeled this approach a ‘fitted-parameter analytical model’ [Refs. 2,3: pp. 
45,73-88]. 

The second method for developing the attrition-rate coefficients has been 
labeled an ‘independant analytical model’ by S. Bonder [Refs. 2,3: pp. 47,73-88]. The 
concept for developing the coefficients is to consider a single firer shooting at a single 
target (which is stationary and does not fire back). These coefficients are developed 
under perfect conditions with this method. There has been a large amount of work 
done in the development of acceptable attrition-rate coefficients. Taylor provides a 
good list of references in his book Force-on-Force Attrition Modeling [Ref. 2]. 

5. Modern Conditions of Warfare 

In 1914 when F. W. Lanchester developed his original models, he was 
attempting to model the combat attrition process so that he could learn more about 
the dynamics of combat under ‘modern conditions’. The phrase modern conditions 
Was an important one then as it is today. In earlier sections Lanchesters models were 
introduced and the assumptions that applied to them were given. The assumption sets 
are what tailor the model to the conditions. Some of these original assumptions no 
longer seem appropriate to model what we would term ‘combat under modern 
conditions’. The next section discusses some of the shortcomings of the original 
Lanchester Models, but before that it is important to lay out a solid definition for 
combat under modern conditions. 

Today the conditions on the combat battlefield are very complex. First, the 
air and land battles are thought of as one battle and given the title “AIR-LAND 
BATTLE’. The ground forces commander from the battalion level on up, has a 
combination of different combat units under his direct command. An Army Division 
Commander has Infantry, Armor, and Artillery Battalions all assigned to his command. 
The battalions which are assigned are not ‘pure’ units, they too have a combination of 
combat units assigned. This concept has been adopted by the Army as the way to 


fight the next war. The Army has organized under the ‘combined arms’ concept and 
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believe that placing a mix of supporting combat forces together on the battlefield will 
be a force multiplier. 

The Navy has developed much the same type of conditions at sea. The 
Carrier Battle Groups (CBG) are combinations of many types of naval combat power 
and therefore brings a variety of forces together to fight a battle at one point under the 
Officer-in-Tactical-Command (OTC). 

6. Shortcomings With The Original Lanchester Models 

Lanchester’s original models provided a very simplified model of modern 
warfare which has served as the basis for a great deal of further work in this area. 
There are many shortcomings of the original models which have been identified and 
have stimulated much of the work done in the years following their development. 
Lanchesters models have been popular because of their simplicity and the ease of 
tracing through the mathematical computations. In this section I will discuss some of 
the shortcoming which have been identified in applying Lanchesters original laws to 
model today’s combat. 

First, not all of today’s combat scenarios can be neatly placed into the 
‘square’ or ‘linear’ attrition models. The most often referenced and easily explained 
example is a guerilla ambush scenario. In this case it is assumed that the guerilla force 
is heavily camouflaged and the force being attacked 1s in the open. It is apparent that 
the attacked force would suffer a ‘squared’ attrition rate from the aimed fire of the 
guerilla force which has clear line of sight of the entire force under attack. On the 
other hand, the defending force would not have clear line of sight of the attacking 
guerilla force and therefore would not be properly modeled by a ‘square’ attrition 
model, but rather the ‘linear’ attrition model. However, even this does not properly 
model the scenario. With time the defending force will be able to locate the attacking 
force well enough that their fire should be considered as ‘aimed fire’ and modeled by a 
‘squared’ attrition model. This illustrates the dynamic conditions on the battlefield 
which must be considered. [Ref. 4] 

One of the main assumptions of Lanchester’s laws is that the two opposing 
forces are comprised of homogeneous units (i.e., Infantry vs.. Infantry). This is rarely 
the case on the battlefield today. Today the doctrine for a land battle is the Air-Land 
Battle Doctrine, and is written for combined arms units. The combined arms concept 
was introduced earlier, but warrants reemphasis. The force that must be modeled 


today combines several different combat units under one commander to fight a battle. 


An additional consideration is the close relationship between the ground force and the 
close air support provided by the Air Force which adds another unit type into the 
model. The same type of unit composition is currently found in the Navy Carrier 
Battle Groups. The need for a model which allows for several different types of units 
(the heterogeneous force model) 1s obvious. The use of heterogeneous forces requires 
the model to consider allocating fire power on a variety of targets. The factors that 
enter in the allocation process are; the attrition rate of each.type of weapon for each 
of the opposing forces weapons types, the number and type of enemy targets, and the 
relative significance of each of the enemy targets. 

Another assumption of the original Lanchester theory is that all of the forces 
of both sides are committed to battle at the beginning and there are no reinforcements 
available. Furthermore, there is no allowance for the possibility of withdrawing some 
or all of a unit’s forces. These factors are definitly consider by the commander on the 
battlefield today. 

Separation of forces on the battlefield is ignored in the original models. They 
assume that the units remain fixed in location and are always within range of the 
opposing forces weapons. To better model the combat process, a movement factor 
which would affect the weapons’ ranges and the accuracy of each weapon on the 
battlefield would be required. 

The original model is designed to model a ‘fight to the finish’ type of combat. 
This definition for battle termination, of fighting to the last man, is not realistic on 
today’s battlefield. There is some point in a battle where the defender will decide to 
withdraw or surrender before the entire force is destroyed. [Ref: 4] 

With this introduction to Lanchester Equation type models as a basis, the 


next chapter will present the two models used for this research. 
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Il. THE MODELS 


A. GENERAL 

Two models were used for this research, one is a deterministic model and the 
Other is a stochastic computer simulation model. The deterministic model used was 
developed at the Naval Postgraduate School by Professor Paul H. Moose and 
Professor John M. Wozencraft. Their model is a Generalized Lanchester Equation 
type model. It will be referenced throughout the remainder of this thesis as the M/W 
model or the Moose/ Wozencraft model. 

The stochastic computer simulation model which was utilized for this research is 
a Monte Carlo Simulation written by Professor Don E. Harrison, Jr. at the Naval 
Postgraduate School. It will be referenced throughout the remainder of the thesis as 
the C/S model or the computer simulation model. 

In the following sections, each of these models will be presented to provide a 


general understanding of there characteristics. 


B. MOOSE/WOZENCRAFT MODEL 

Paul H. Moose is a professor in the Electrical Engineering Department at the 
Naval Postgraduate School and John M. Wozencraft was a professor with the 
Electrical Engineering Department prior to his retirement. Both have been associated 
with the academic group for the Joint Command, Control, and Communications 
Curriculum, and have been interested in the decision and control problems a modern 
day military commander faces. They wanted to develop a better understanding of how 
command and control decisions affect battle outcome. To do this, an analytical model 
of the attrition process that would provide an adequate representation of the modern 
conditions of combat was needed. Such a model required modeling a variety of forces 
as an aggregate force fighting on the battlefield. [Ref. 5] 

As presented in chapter one, the military commander in today’s environment has 
several types of units available to him to fight a battle. The question of optimum 
resource allocation involves the use of a variety of these units. For the model to be 
useful it had to be easy to interpret and understand. Lanchester Equation type models 
were easy to interpret, but were limited in their ability to model a variety of forces. 
This motivated their research into the dynamics of a Generalized System of Lanchester 


Equations of the type shown here in Equation Set 2.1: 
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x(t) = E (X:,y,t): WHETC es ie 
(eqn 2.1) 
*i(t) = Gi(Xy;,t): where j= 1,2,........M, 


where N components of a non-homogeneous X-force engage M components of a Y- 
force. As shown in Equation Set 2.1, the number of combatants of type x; is a 
function of their own forces, the total Y force (the sum of the Yj forces) available to 
fight, and time. Similarly, the rate of change of a Y; force is a function of it’s own 
force level, the total X force (the sum of all x; forces) available and time. 


The following set of equations represent the Moose/Wozencraft model: 


X(t) = -ux: (tye ai i(t)y;(t)- 2.05 jy) + tj] where i=1,2,.....,.N, 
(eqn 2.2) 
yj(t)- “Liciix i(ty;(t)- 24) x(t) + s; where j= 1,2,....,M, 


where 
1) The x; and y; Tepresent Blue and Orange forces of different types 
2) The u; and v, are self attrition coefficients 


J 
3) The aij and Cj; are area fire coefficients 
4) The bi; and di are aimed fire coefficients 
5) Ther; and sj are resupply coefficients 


6) Xi; 2 QO for alli andj 


7) 2 X=X%, » Vie 

The M/W model is designed as an NxM, heterogeneous force, combat attrition 
model. The significance of this fact is that it allows for closer modeling of the modern 
combat force discussed in Chapter One. 

In the M/W model, Equation Set 2.2, several things are different from the 
original Lanchester Equation models presented in Chapter One. One major difference 
is that the M/W model allows for modeling N type X forces and M type Y forces. The 
Moose/Wozencraft model treats the attrition of a given force as one process where 
Lanchester’s original work broke the attrition process into two models, aimed and area 
fire models. With that fact in mind then, the rate of change of a force x; is a function 


of four separate factors; self-attrition, area fire attrition, aimed fire attrition and 
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resupply. In the case of determining the attrition to force x; , the self-attrition is 
modeled using a self-attrition coefficient -u; which is multiplied by the current x; force 
level. This yields the attrition to force x; which is attributed to ‘self losses’ such as 
disease and defections. The area fire attrition 1s modeled similarly to Lanchesters 
original model, with the added capability to account for all the possible combinations 
of Yj forces which can cause casualties to x; by means of area fire. Thus an attrition- 
rate coefficient for each possible combination of forces engaging in battle must be 


available. Each Yj force has an area fire coefficients -a;. associated with it. The sum of 


i 
the attrition-rate coefficients (aj;) multiplied by the Yj force level can then be thought 
of as an aggregate area fire attrition rate which multiplies the x; force level. The aimed 
fire is very similar. It is also similar to the original Lanchester model with the added 
capabilities to model NxM force battles. The aimed fire is modeled by multipling each 
Yj force by an attrition-rate coefficient “bj and summing them up. This is then the 
aggregate aimed fire attrition rate. Finally the last factor modeled is the resupply of 


each x; force with a resupply variable r;. The sum of the self-losses, area fire losses, 


i: 
aimed fire losses and resupply is equal to the change in the x; force. 

The model determines the change in a Yj force in exactly the same manner. The 
attrition rate coefficients are different of course, but the process is the same. The self- 
attrition coefficient for Y; 1S “Vje The area-fire attrition coefficients are given as Ci and 


the aimed fire coefficients are given as -d;;, The resupply variable for Yj is $ 


The relative usefulness of this mala for modeling the combat in process 
on today’s battlefield is apparent from this discussion. First, it provides the capability 
to model a battle involving two heterogeneous forces which we saw is necessary to 
analyze the optimum force allocation issue. Secondly, it is a deterministic model, so it 
is easy to understand and interpret. It has many similarities to the original Lanchester 
Equation models which are widely studied and understood. 

The results of Professor Moose and Professor Wozencraft’s research into the 
dynamical properties of this model are provided in a paper which they are submitting 
to the Military Operations Research Society for publication titled Characteristic 


Trajectories of Generalized Lanchester Equations [Ref. 6]. 


C. STOCHASTIC COMPUTER SIMULATION MODEL 
The stochastic model is a Monte Carlo Simulation designed to test the 


Moose/Wozencraft model. The basis of this model was designed by Don E. Harrison, 
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Jr., a professor at the Naval Postgraduate School. Professor Harrison was asked to 
design a stochastic model of the attrition process which included the same factors as 
the Moose/Wozencraft model. The combination of his efforts and some additional 
code to run the program for a given number of replications and calculate statistical 
data resulted in the model shown in Appendix B entitled Computer Simulation model 
(GIS). 
The C/S model will be presented in the following sections: 
1) Input parameters 
2) The Combat Cycle 
3) The Output 
The material in each section is presented in a condensed users manual form to provide 
the reader with an understanding of the model and the capability to exercise the model 
if desired. 
1. Input parameters 

The computer code is written in FORTRAN/77 and uses an exec file to define 
all the input files, output files, load and execute the program. 

The input file contains all the required input data. Figure 2.1 shows the 
contents and format of the input data file. The first line of the file is a heading which 
is read by the program and allows the user to label and keep track of which test data 
set is being utilized. Next the data required for the x; forces is given. On the second 
line the number of types of X forces is given (the variable NI YPX in the program), 
this corresponds to the N value in the M/W model which specifies the number of types 
of X forces. Then the initial force level for each X; force is given, their locations, their 
corresponding self-attrition coefficients (u;) and the resupply variables (r;). Next, the 
attrition-rate coefficients are provided in matrix form. The aimed fire attrition-rate 
coefficients ‘d:.’ are given in lines three through five. Area fire attrition-rate 


1 

coefficients ci are given in lines six through eight. Then the required information for 
the Y; forces is listed. Once again, the number of type of Y; forces is given (the 
NTYPY variable in the program) which corresponds to M in the M/W model. Then 
each Y; initial force level, location, self attrition-rate coefficients Vy and resupply 
variables ‘Sf are given. Lines ten through twelve contain the aimed fire coefficients 
bi, and lines thirteen through fifteen contain the area fire coefficients ai Cite 
sixteen of the file contains seed numbers for the pseudo-random number generators. 
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01/16/87 TEST OF (2X2), 1 UNSTABLE ROOT $ 
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Figure 2.1 Input Data File. 


The program is designed to be interactive before it begins the combat cycle 
which allows the user to change the heading of the input file if so desired. This also 
enables the user to ensure that the input file is the one desired. 

2. The Combat Cycle 

The combat cycle is a stochastic process. During the combat cycle the 
combatants are chosen at random and the probability of self-loss, aimed fire loss, and 
area fire losses are all tested by use of a pseudo-random number between 0.0 and 1.0. 
The cycle begins by randomly selecting which force, X or Y, will fire. Once the force 
has been selected then a combatant from the force is selected. For the purpose of 
presentation, say that x,(1) (that is combatant one of force type x1) is selected to be 
the firer, then the first thing that occurs in the cycle is a check to see if combatant 


x,(1) is a self loss. This is accomplished by use of a random number (between 0.0 and 
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1.0) which is compared against the appropriate self attrition coefficient uy. Tf ihe 
combatant 1s not a self loss then a Y force combatant is selected at random as a target, 
say y>(2) (the second combatant of force type V2) as an example. Once the firer and 
target have been identified, the next step tests the effectiveness of the ‘shot’ by using a 
random number to test whether the target is killed. This process involves first 
determining if the target is killed, which is done by testing whether a random number is 
less than the total kill probablity (which is the sum of the aimed fire and area fire 
attrition-rate coefficients, or their kill probabilities). If the random number is less, then 
the target has been killed and if not the shot missed. If the target was killed then this 
same process 1s followed to determine whether it was a loss due to aimed fire or area 
fire. Required data is tabulated on the number of combatants of type Yj Which are 
killed by type x; forces and each force is resupplied. The combat cycle is repeated in 
this manner until all combatants have had an opportunity to fire during each timestep. 
This process 1s continued for each timestep until one of two conditions occurs: 

1} One force reaches a break point which has been defined as NXSTOP and 
NYSTOP. This simulates the level of attrition a commander will suffer before 
pulling back or surrendering rather than fighting to the last man. 

2) The combat cycle has been repeated for a given number of timesteps, which 1s 
specified in the input data file. 

3. The Output 

The C/S model was designed to provide several output files for data analysis. 
The model generates two files which provide a very detailed timestep-by-timestep 
recording of the combat cycle outcome. These files require a great deal of storage 
space and therefore during the interactive portion of the program the user is asked if 
they want these files printed in their complete form. It is not recommended that they 
be printed if the user desires multiple replications because the storage space required 1s 
too large. One of these files provides a detailed listing of each timestep results 
including a summary of the force levels, the number of losses by each type of fire and 
the number of combatants that were resupplied. An example of this file’s output is 
shown in Figure 2.2. The other file contains similar listings with the number of dead 
combatants of each type provided as well. An example of this file’s output is shown in 
Figure 2.3. 

The remainder of the files are designed to provide the force levels of each force 


(1.€., X; and y,) for each timestep. This provides the data which can be used to plot 
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02/02/87 TEST OF (3X3), 
# X TYPES NUMBER EACH LOCATIONS SELF-LOSS RE-SUPPLY 

3% 200 150 175 1.0 1.0 1.0 0.005 0.010 0.020 a 4 ll 
KILL PROBABILITIES FOR FORCE X(I) SHOOTS Y(J) 
PK( 1,1) PK(1,2) PK(1,3) PK(2,1) PK(2,2) PK(2,3) PK(3,1) PK(3,2) PK(3,3) 
0.06000 0.02000 0.02000 0.00000 0.00000 0.00000 0.03000 0.01000 0.02000 
AREA FIRE COEFFICIENTS FORCE X(I) SHOOTS Y(J) 


CC(1,1) cC(1,2) CC(1,3) CC(2,1) CC(2,2) CC(2,3) CC(3,1) CC(3,2) CC(3,3) 
0.00000 0.00000 0.00000 0.00000 0.00030 0.00020 0.00000 0.00000 0.00010 


# Y TYPES NUMBER EACH LOCATIONS SELF-LOSS RE-SUPPLY 
3 150 125 225 1.0 1.0 1.0 0.003 9.020 0.030 ‘+ 5 14 


KILL PROBABILITIES FOR FORCE Y(I) SHOOTS X(J) 
PK(1,1) PK(1,2) PK(1,3) PK(2,1) PK(2,2) PK(2,3) PK(3,1) PK(3,2) PK(3,3) 
0.05000 0.01000 0.05000 0.00000 0.00000 0.00000 0.02000 0.01000 0.03000 
AREA FIRE COEFFICIENTS FORCE Y(I) SHOOTS X(J) 


AA( 1,1) AAC 1,2) AAC1,3) AAC2,1) AA(2,2) AAC 2,3) AAC3,1) AAC3,2) AA(3>3) 
0.00000 0.00000 0.00000 0.00000 0.00020 0.00030 0.00000 0.00000 0.00010 


DOUBLE PRECISION RANDOM VARIANT SEEDS INTEGER SEED VALUES 
0.34521579D+09 0.45635761D+08 0.89342761D+08 1209308051 1194034571 1196766258 


# NTXS NTX NX1l NX2 NX3 RKX1 RX2 RX3 NTYS NTY NYL NY2 NY3 RY1 RY2 RY3 
0 525 525 200 150 175 @ 4@ 11 500 500 150 125 225 %@ 5 14 
TIMESTEP 1 
AIMED FIRE KILLS OF Y"S BY X FORCE 


D(1,1) O(1,2) OD(1,3) O2,1) Dl2,2) Dl2,3) O(3,1) O(3,2) O0(3,3) 
3 0 2 0 0 0 1 iL 0 


AIMED FIRE KILLS OF X"S BY Y FORCE 


B(l,1l) Bl,2) Bll,3) Bl2,1) Bl2,2) Bl2,3) BI3,1) BI3,2) B3,3) 
2 1 1 0 0 0 2 ) 2 


AREA FIRE KILLS OF Y"S BY X FORCE 


C(l,1) Cll,2) Cll,3) Cl2,1) Cl2,2) Cl2,3) Cl3,1) Cl3,2) Cl3,3) 
0 0 0 0 5 7 0 0 2 


AREA FIRE KILLS OF X"S BY Y FORCE 


A(1l,1) Al1,2) Al1,3) Al2,1) Al2,2) Al2,3) Al3,1) Al3,2) Al3,3) 
0 0 0 0 Z 5 0 0 3 


# NTXS NTX NX1 NX2 NX3S DX1 OX2 DX3 NTYS NTY NY1 NY2 NY3 OY1 DY2 DY3 
1 525 525 200 150 175 5 5 13 500 500 150 125 225 4 9 20 


Figure 2.2. Example File 6 Output. 


2/ 


TOTAL X FORCE LOSS WAS 23, RESUPPLY WAS: 
AIMED FIRE LOSSES: 4 1 3 
AREA FIRE LOSSES: 0 2 8 
SELF LOSSES: 1 2 2 
TOTAL Y FORCE LOSS WAS 33, RESUPPLY WAS: 
AIMED FIRE LOSSES: 4 1 2 
AREA FIRE LOSSES: 0 5 9 
SELF LOSSES: 0 3 9 
1 525 521 199 149 173 4 @ 11 500 490 150 121 219 4 5 14 


AVERAGE ATTRITION COEFFICIENTS FOR Y FORCE 


AA(1,1) AA(1,2) AA(1,3) AA(C2Z,1) AA(2,2) AA(2,3) AA(3,1) AA(3,2) AA3,3) 
0.00000 0.00000 0.00000 0.00000 0.00011 0.00023 0.00000 0.00000 0.00008 


AVERAGE ATTRITION COEFFICIENTS FOR Y FORCE 


BB(1,1) BB(1,2) BB(1,3) BB(2,1) BB(2,2) BB(2,3) BB(3,1) BB(3,2) BBI3,3) 
0.03500 0.02333 0.02000 0.00000 0.00000 0.00000 0.02333 0.00000 0.02667 


AVERAGE ATTRITION COEFFICIENTS FOR & FORCE 


CC(1,1) CC(1,2) CC(1,3) CC(2,1) CC(2,2) CC(2,3) CC(3,1) CC(3,2) CC(3,3) 
0.00000 0.00000 0.00000 0.00000 0.00027 0.00021 0.00000 0.00000 0.00005 


AVERAGE ATTRITION COEFFICIENTS FOR X FORCE 


DD(1,1) DDO(1,2) DD(1,3) DD(2,1) DD(2,2) DD(2,3) DD(3,1) DD(3,2) DDI3,3) 
0.05000 0.00000 0.02222 0.00000 0.00000 0.00000 0.01905 0.02286 0.00000 


STOPPED BY PROGRAM AT NSTOP = 1 





Figure 2.2. Example File 6 Output. (cont’d.) 


force level trajectories, which allows for ease in performing a comparison of the two 
model's results for a given input set. 

Now that each of the models has been presented, chapter three will discuss the 
design of the experiment, and provide samples of the output files used for plotting the 


force level trajectories. 
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02/02/87 TEST OF (3X3), 
# NTXS NTX NX1 NX2 NX3 DOX1 OX2 DX3 NTYS NITY NY1 NY2 NY DY1 OY2 DYS 
0 525 525 200 150 175 0 0 0 500 500 150 125 225 0 0 0 


D,B,C,A BY ROWS IN DESCENDING ORDER 
3 0 2 0 
2 1 1 0 
0 0 0 0 
0 0 0 0 
TOTAL X FORCE LOSS WAS 23, RESUPPLY WAS: G¢ | 
AIMED FIRE LOSSES: 4 1 3 
AREA FIRE LOSSES: 0 2 8 
SELF LOSSES: 1 2 2 
TOTAL Y FORCE LOSS WAS 33, RESUPPLY WAS: ¢ 5 14 
AIMED FIRE LOSSES: o 1 2 
AREA FIRE LOSSES: 0 5 9 
Seun eOsseo. 0 3 Ms 
L 525..521..199 149 173 5 5 13 500 490 150 121 219 «© 2 


1 
1 
1 
1 


A,B»,C,D BY ROWS IN DESCENDING ORDER 
AV 0.00000 0.00000 0.00000 0.00000 0.00011 0.00023 0.00000 0.00000 0.00008 
AV 0.03500 0.02333 0.02000 0.00000 0.00000 0.00000 0.02333 0.00000 0.02667 
AV 0.00000 0.00000 0.00000 0.00000 0.00027 0.00021 0.00000 0.00000 0.00005 
AV 0.05000 0.00000 0.02222 0.00000 0.00000 0.00000 0.01905 0.02286 0.00000 
STOPPED BY PROGRAM AT NSTOP 1 





Figure 2.3 Example File 7 Output. 
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It. RESEARCH METHODOLOGY 


A. OBJECTIVE 

The primary objective of this research is to determine whether a deterministic 
combat attrition model 1s a good approximation for the mean result of a stochastic 
niodel of the combat attrition process. With that objective in mind, an experiment was 
designed to preform a comparison of two such models. The models used in this 
experiment are the Moose/Wozencraft model (M/W model), which is a deterministic 
model, and a Monte Carlo Simulation (C/S model) which were introduced in chapter 
two. A secondary objective is to see if the two models behave with similar 
characteristics as the results of Professor Moose and Professor Wozencraft’s research 
[Ref. 6]. 


B. EXPERIMENTAL DESIGN 
The experiment is designed to run-the C/S model for ninety-nine timesteps over 
thirty replications, to generate force level data and calculate the resulting average 
attrition rate coefficients. These average coefficients are then used as the input 
parameters to the M/W model which then provides the calculated force levels per 
timestep. This force level data is used for analysis. The method for analysis which has 
been chosen 1s to plot the resulting force level data over time to provide for an easier 
comparison than sorting through the enormous amount of data generated for each test 
case. | 
This experiment was conducted for cases involving 1x1, 1x2, 2x2, and 3x3 force 
level scenarios. Examples of the data collected and the resulting force level trajectory 
plots are provided in the next section. 
1. Model Verification 
Before beginning the experiment using these two models, each of the models 
was exercised and tested heavily to verify that they ran properly and provided the 
desired output as they were designed. A brief description of the verification process for 
each model will be given. 
a. Computer Simulation Model (C/S model) 
The first step in verifying the C/S model involved a procedure of several 


runs using different input data sets to test the program to ensure that: 
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1) The input data and program were compatible and the input variables are 
properly read. 

2) All the computational code was performing properly (1.e., the average attrition 
rate coefficients and the average force levels). 

3) That each of the test conditions was functioning as designed (1e., the 
condition for ending at the specified force levels XSTOP and YSTOP). 

4) That the output was presented as desired for the step-by-step analysis and the 
plotting data. 

5) Test the resulting average attrition-rate coefficients with the input probabilities 
of kill which correspond to these coefficients. 

During this portion of the research it was interesting to note that the C/S 
model was designed to allow each shooter one shot per time interval, but the shooter 
could select a target to fire at which had already been killed during that time interval. 
This fact seemed to be a fairly realistic condition that could occur on the battle field 
where two combatants fire at the same target instantaneously or within a fraction of a 
second of each other. Therefore, this condition was left in the model. 

b. Moose] Wozencraft Model (M/W model) 

A computer program was designed to run the M/W model to allow for an 
easy comparison of the two models. By utilizing a computer program it allowed for 
the same timestep by timestep analysis as the C/S model. 

After designing the FORTRAN program, shown in Appendix D, to run the 
Moose/Wozencraft model (M/W model) it was tested to verify that the following 
conditions were met: 

1) The input data and the program were compatible and all variables were read 
properly. 

2) All computational code was performing as designed. 

3) The resulting output files provided the desired data for analysis. Since the 
M/W model is a deterministic model, the same output will be obtained each 
time for the same set of input data parameters. 

Each of the conditions listed above was verified for simple input data sets for each of 
timencases. 
2. Data Generation Procedure 
The procedure followed in conducting this experiment was a three step process 
for each unique case tested. The first step involved running the C/S model for ninety- 


nine timesteps over thirty replications which would generate the following data files: 
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1) The calculated average attrition-rate coefficients which are later used as input 


to the M/W model. An example file is shown in Figure 3.1. 


01/27/87 TEST OF (1X1) 
AVERAGE ATTRITION COEFFICIENTS FOR Y FORCE 


AA( 1,1) AAC1,2) AAC1,3) AA(Z,1) AAC2Z,2) AA2Z,3) AAC3)1) 1 AACS, 2) AA(3,3) 
0.00012 0.00000 0.00000 0.00001 0.00000 0.00000 0.00000 0.00000 0.00000 


AVERAGE ATTRITION COEFFICIENTS FOR Y FORCE 


BB( 1,1) BB(1,2) BB(1,3) BB(2,1) BB(2,2) BB(2,3) BBI(3,1) BB(3,2) BBI3,3) 
0.04857 0.00000 0.00000 0.00142 0.00000 0.00000 0.00000 0.00000 0.00000 


AVERAGE ATTRITION COEFFICIENTS FOR X FORCE 


CC(1,1) CC(1,2) CC(1,3) CC(2,1) CC(2,2) CC(2,3) CC(3,1) CC(3,2) CC(3,3) 
0.00027 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 


AVERAGE ATTRITION COEFFICIENTS FOR X FORCE 


0O(1,1) DO(1,2) DDO(1,3) DD(2,1) DDO(2,2) DDC 2,3) OD(3,1) DO(3,2) DD(3,3) 
0.08648 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 


STOPPED BY PROGRAM AT NSTOP = 





Figure 3.1 Example Average Attrition-rate Coefficients File For C/S Model. 


2) A file containing the average force level for each force, the average force level 
plus one standard deviation, and the average force level minus one standard 
deviation is provided for the aggregate force levels X (2.x) and Y DY; -) per 
timestep as well as for each of the individual units x; and Yj: An ail file 
is shown in Figure 3.2. 

3) The average force level and standard deviation of each force level per time 
step. An example of these results is shown in Figure 3.3. 

The second step involves using the average attrition-rate coefficients calculated 
by the C/S model as input parameters to the M/W model. The output generated by 
the M/W model for each set of input data is the following: 

1) The force level for each of the forces x; and y; as well as the total force levels 
x (1X) and Y (Lj) for each time step. An example of this output file 1s 


shown in Figure 3.4. 


ar 


01/16/87 TEST OF (2X2), 1 UNSTABLE ROOT 





Figure 3.3 Example File of Force Level and Standard Deviation For C/S Model. 


The third step involves plotting the force level data generated by each of the 
models to perform the comparison. The plots which were utilized for the comparison 
are listed below with examples provided: 

1) <A plot of the total X force vs. Total Y force level trajectories for both models 


is done first. This plot provides for a macro level comparison of the resulting 


aS 


01/27/87 TEST OF (1x2), CASE 4 


ae Xy a 

Q 200.00 300.00 
1 196.29 299.19 
2 192.88 298.75 
3 189.73 298.61 
4 186.82 298.76 
5 184.12 299.13 


Figure 3.4 Example xy> Force Level Data File For M/W Model. 


force levels over time (1.e., the Force Commanders level). This plot is two 
dimensional, plotting the X force level vs.. Y force level. The time variable 1s 
shown by marking the force level every fifth timestep. This plot will be used 
to analyze the two aggregate force models results. An example of the total 
force level trajectory plot is shown in Figure 3.5. 

2) Plots of the x; force level vs.. time and the Yj force levels vs. time for both the 
C/S and the M/W models are given. This provides a comparison of the battle 
outcome for the two models at the micro level. An example of the x; (or yj) 
plot is shown in Figure 3.6. 

The plots introduced in this chapter are used extensively for conducting the 
comparison of the two models. Chapter Four will present the results of the experiment 


and discuss the comparisons. 
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Figure 3.5 Example Total Force Level Trajectory Plot. 
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Figure 3.6 Example plot of x; force level trajectories. 
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IV. EXPERIMENTAL RESULTS 


A. INTRODUCTION 

The primary objective of this experiment was to test whether the deterministic 
Moose; Wozencraft model was a good approximation (or prediction) for the mean 
result of the stochastic computer simulation. The secondary objective for this 
experiment was to test whether the two models’ force level trajectories behaved in the 
expected manner described by Professors Moose and Wozencraft. The results of ten 
cases are presented graphically showing the two models’ results with a bracket of plus 
and minus one standard deviation as a measure of uncertainty. Each case will be 
discussed in the context of the primary and secondary objectives in this chapter and 


their results shown graphically in Appendix E. 


B. THE FOUR MODES OF BEHAVIOR 

Before proceeding with the presentation of the results, it 1s important to briefly 
explain the four trajectory behavior modes which will be referenced in each of the 
cases. Professors Moose and Wozencraft found that the force level trajectories 
displayed one of four distinct modes of behavior. A mode of behavior was the result of 
the existence, or nonexistence, of equilibrium points in the state space (positive first 
quadrant). The four modes related to the existence of one stable equilibrium, one 
unstable equilibrium, two equilibria (where one is stable and one is unstable), or no 
equilibria. These four modes of behavior are illustrated in Figure 4.1. [Ref. 6] 

The expected behavior for any force level trajectory when one stable equilibrium 
point exists in the state space is shown in Figure 4.1 (a). The trajectory will be 
attracted to the stable equilibrium and remain there, since dx/dt=dy/dt=0 at the 
stable equilibrium. 

The expected behavior for any force level trajectory when one _ unstable 
equilibrium point exists in the state space is shown in Figure 4.1 (b). As the 
trajectories near the equilibrium they are repelled or turned away from the equilibrium 
point. The presence of an unstable equilibrium creates a division of the state space 
into two regions. The dividing line was called a separatrix. In scenarios which exhibit 
this mode of behavior the outcome depends on which side of the separatrix the initial 


force levels begin. 
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Separatrix 





Figure 4.1 The Four Trajectory Modes Of Behavior. 
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The expected behavior for any force level trajectory when two equilibria exist 
(one stable and one unstable) 1s shown in Figure 4.1 (c). Once again the unstable 
equilibrium creates a separatrix dividing the state space into regions. The stable 
equilibrium will attract trajectories as discussed earlier. This particular mode of 
behavior has two possible outcomes for a battle, one of the forces can win or they will 
fight to a draw. The outcome depends on the initial force levels for each force. 

The expected behavior for any force level trajectory when no equilibrium point 
exists is shown in Figure 4.1 (d). The trajectories will approach, and eventually cross 


one of the axes, which results in a force being eliminated. 


C. DISCUSSION 

This numerical experiment was designed to test a combat scenario involving 
different opposing-force composition. (N, x M,,) cases involving N, types of X forces 
VS. My types of Y forces are tested in this experiment. An example is a 1x2 case which 
involves one X type force vs. two Y type forces. The input data sets were chosen to 
test the two models for the four distinct modes of behavior that Professors Moose and 
Wozencraft identified [Ref. 6]. The computed results are presented in sections which 
are dedicated to a specific (N x M) scenario. 

Appendix E contains the input data sets used for each model, and the resulting 
force level trajectory plots for each of the cases. In each case the input data for the 
C/S model was developed to test one of the modes of behavior and the M/W model 
input data were the calculated attrition rate coefficients calculated by the C/S model. 
Several trajectory plots were done for each of the cases. There is a plot of the total 
aggregate force level trajectory which compares the total X force vs. the total Y force. 
This aggregate trajectory shows the battle outcome. An example of this plot is shown 
for case six in Figure 4.2. 

The result of the computer simulation model, which is the mean of thirty 
replications, is always shown as a solid line with squares marking the force level every 
fifth time step. The Moose/Wozencraft model results are shown as a dashed line with 
circles marking the force level every fifth time step. Additionally, two solid lines are 
shown with no markings. These are the curves which show plus and minus one 
standard deviation from the C/S model results. They serve as a measure of uncertainty 
for the comparison. The aggregate force level plot provides a comparison of the two 


models total force level results to determine whether they do have the same general 
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Figure 4.2 Total Force Level Trajectory For Case Eight. 
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behavior and whether the deterministic model 1s a good approximation for the result of 
the stochastic model. Each of the aggregate forces (1e., X = x) are plotted as a 
function of time to provide a clearer comparison of the two models results as the battle 
progresses. An example of these plots are shown in Figures 4.3 and 4.4. This provides 
another comparison between the two models’ results for each aggregate force. In 
addition to the aggregate force level plots, the results for each of the separate units are 
plotted as a function of time. Examples of these plots are provided in Figures 4.5 
through 4.8. This provides a comparison of the two models at the individual unit level. 

The method of comparison used for this analysis was a graphical presentation of 
the force level trajectories rather than an extensive numerical analysis of the resulting 
data. Each case will be presented in a similar manner. The cases will be briefly 
described and then a list of observations given. Specific comments are made about 
whether the M/W trajectory was close to the C/S trajectory and whether the expected 


behavior was observed. 


D. (1X 1) CASE COMPARISONS 

The experiment was conducted for cases one through four using one X force and 
one Y force, which are labeled as 1xl cases. The results of each of these cases are 
shown in Figures E.1 through E.12 in Appendix E. The force level trajectories for the 
M/W model were very close to the mean force level trajectories for the C/S model in all 
four cases. 

1. Case One 

Case one was designed to test the two models involving a 1x1 scenario for the 

mode of behavior where no equilibrium exists. The following observations are noted 
from the results shown in Figures E.1 through E.3: 

a) The resulting force level trajectories for the M/W model and the C/S model 
were very similar, both displaying the same trajectory shape. 

b) The M/W model trajectory was well within the one standard deviation 
boundaries, for the complete battle. In fact, the two curves over lapped 
during portions of the battle. This case resulted in the X force winning the 
battle. 

2. Case Two 
Case two was designed to test the mode when one stable equilibrium exists in 
the vicinity of (200,200). The following observations are noted from the results shown 
in Figures E.4 through E.6: 
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Figure 4.3 X Force Level Trajectory Over Time For Case Eight. 
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Figure 4.4 Y Force Level Trajectory Over Time For Case Eight. 
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Figure 4.5 X1 Force Level Trajectory Over Time For Case Eight. 
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Figure 4.6 X2 Force Level Trajectory Over Time For Case Eight. 
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Figure 4.7 YI Force Level Trajectory Over Time For Case Eight. 
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Figure 4.8 Y2 Force Level Trajectory Over Time For Case Eight. 


47 


a) 


b) 


Cc) 


The M/W trajectory and the C/S trajectory were very similar and they display 
the same trajectory shape. 

The M/W model trajectory was outside the established boundaries for the 
initial portion of the battle. It is noted, however, that the area bounded by 
one standard deviation is extremely small in this region. The M/W trajectory 
entered the bounded area as the two curves turned. Both models reached an 
equilibrium point in the vicinity of (200,200) where the battle would end in a 
draw. 

Examination of the individual unit force level trajectories in Figures E.5 and 
E.6 shows that the M/W trajectory was within the established bounds for each 


force as a function of time. 


. Case Three 


Case three was designed to test the 1xl scenario for the mode of behavior 


where one unstable equilibrium exists in the vicinity of (200,200). The following 


observations are noted from the results shown in Figures E.7 through E.9: 


a) 


b) 


The force level trajectories of both models displayed the expected behavior for 
a case where one unstable equilibrium exists. Each of the models trajectories 
approached the equilibrium and then turn away from it. 

Once again the M/W trajectory was well within one standard deviation of the 


C/S trajectory. 


. Case Four 


Case four was designed to test a Ixl scenario where the mode of behavior 


involves two equilibria in the vicinity of (200,200)(U) and (366,100)(S). The following 


observations are noted from the results of this case shown in Figures E.10 through 


Ear? 
a) 


b) 


E: 


The results show the trajectories of both models display the expected behavior 
approaching the unstable equilibrium and then turning away. 
The M/W model was always well within the one standard deviation 


boundaries for this case. 


(1X2) CASE COMPARISONS 


One case was developed to test an unbalanced force level scenario involving one 


type X force and two type Y force. The purpose of this case was to ensure that both 


models perform properly for an unbalanced number of forces engaging in the battle. A 


1x2 scenario was chosen for this case and the results are presented next. 
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1. Case Five 


Case five was the one unbalanced scenario developed and the mode of 


behavior tested involved one unstable equilibrium in the vicinity of (300,200). The 


following observations were noted from the results shown in Figures E.13 through 


E.18: 
a) 


b) 


Cc) 


d) 


P 


Both models displayed the expected behavior for an unstable equilibrium. The 
trajectories for both models approached the equilibrium point and then were 
turned away. 

The M/W trajectory was well within the one standard deviation boundaries, 
and had the same shape as the C/S trajectory. It is interesting to note the 
relative speed at which the two trajectories moved away from the unstable 
equilibrium point. The C/S model moved away much faster than the M/W 
model. 

The plot of each force vs. time, shown in Figures E.14 through E.18, show the 
two models results separating over time, and the M/W trajectory leaving the 
bounded area toward the end of the battle. 

This case involving a 1x2 scenario with an unstable equilibrium is similar to 


case three which was a 1x1 scenario involving one unstable equilibrium. 


(2X2) CASE COMPARISONS 


The next step in the experiment was to develop cases which involved two types of 


forces on each side and compare the models results for the four modes of behavior 


discussed. Cases six through nine are 2x2 case scenarios. Their results are shown in 


Figures E.19 through E.46. The results for each of these cases will be discussed 


individually. 


1. Case Six 


Case six is a 2x2 force level scenario designed to test the models when there 


are no equilibria in the state space. The results for case six are shown in Figures E.19 


through E.25. The following observations are noted from the results: 


a) 


The resulting trajectories shown in the Total Force Level Trajectory plot in 
Figure E.19 show both models behaved similarly. In this particular case the Y 
force was the winning force. The expected behavior for no equilibria was 


observed in this case. 
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b) 


Cc) 


d) 


The M/W model results were within the defined one standard deviation 
boundaries during most of the battle, leaving the boundaries only toward the 
end of the engagement. 

Further comparison of the two models can be done using the results plotted 
for the individual forces in Figures E.20 through E.25. From these results it is 
evident that the X force results are the ones where the models differed the 
most. The trajectories all had the same general trends. 

This 2x2 case involving no equilibria is similar to case one which was a Ixl 


case scenario involving no equilibria. 


. Case Seven 


Case seven was designed to test the models results with two type X forces and 


two type Y forces when two equilibria exist. The results of the two models are shown 


in Figures E.26 through E.32. The following observations are noted for case seven: 


a) 


b) 


The resulting trajectories for both models behaved in the expected manner. 
They approached the unstable equilibrium and then turned away and headed 
off. In this particular case the X force is winning, but the results would be 
different if the force levels were such that the initial point was above the 
Separatrix. 

The M/W trajectory had the same shape as the C/S trajectory. The M/W 
trajectory shown in Figure E.26 was outside the established boundaries during 
the initial portion of the battle, but was within the boundaries after they 
turned away from the equilibrium. However, examination of the individual 
units force level trajectories shown in Figures E.27 through E.32 shows the 
M/W model within the boundaries for the X force and just outside the 
boundaries for the Y force. 

This 2x2 case involving two equilibria is similar to case four which was a Ixl 


scenario involving two equilibria. 


. Case Eight 


Case eight was designed to test the two models results for a force scenario 


involving two type X forces and two type Y forces with one stable equilibrium present 


in the state space. The results for case eight are shown in Figure E.33 through E.39. 


The following observations are noted for case eight: 


a) 


The resulting trajectories for the M/W and C/S models both exhibited the 


expected behavior with one stable equilibrium present. They both were 
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b) 


Cc) 


attracted to the stable equilibrium point and remained in its vicinity. The C/S 
model shows more variability due to its stochastic process. 

The M/W trajectory had the same shape and remained within the established 
boundaries for the entire battle. 

This 2x2 case involving one stable equilibrium is similar to case two which was 


a 1x1 scenario involving one stable equilibrium. 


. Case Nine 


Case nine was designed to test the two models results for a 2x2 force level 


scenario when one unstable equilibrium exists. The results for case nine are shown in 


Figures E.40 through E.46. The following observations are noted for case nine: 


a) 


b) 


G. 


The resulting trajectories for each of the models behaved in the expected 
manner for an unstable equilibrium. 

The M/W trajectory had the same general shape, approaching and turning 
away from the equilibrium, but it was not within the designated boundaries 
during most of the battle. Reviewing the results of the individual units shows 
that the two models did not behave similarly as a function of time. 

This 2x2 case involving one unstable equilibrium is similar to case three which 


is a 1x1 scenario involving one unstable equilibrium. 


(3X3) CASE COMPARISONS 


The last force level scenario developed for this experiment involves three type X 


forces and three type Y forces. One case, labeled as a 3x3 case, was developed for this 


experiment. There is no discussion of the two models results for each of the four 


modes of behavior discussed earlier for the 3x3 scenario. Rather, one case was run to 


test the models results for the larger 3x3 scenario. 


1. Case Ten 


Case ten is the one 3x3 case developed for comparing the M/W model results 


with the C/S model results. The following observations are made from the results 
shown in Figures E.47 through E.55: 


a) 


b) 


The M/W trajectory was well within the established boundaries throughout 
the battle and the two trajectories shapes were almost identical. 

The M/W and C/S trajectories for each of the individual units are very simular. 
In several of the cases the two trajectories are overlapping for most of the 
battle. 


>| 


H. SUMMARY OF RESULTS 

The results of the numerical experiment, with ten cases involving Ixl, 1x2, 2x2, 
and 3x3 force level scenarios, have been presented and briefly discussed in this chapter. 
Some general observations from these results will serve as a summary of the 
experiment. 

First, the results for the ten cases support the idea that a deterministic attrition 
model is often a good approximation for the mean result of a stochastic attrition 
model. From the results shown in the previous section, only one of the ten cases 
displayed force level trajectories where the M/W model was significantly different from 
the C/S model. Five of the ten cases showed the M/W trajectories were within the 
established one standard deviation boundary for the entire battle (cases one, three, 
four, eight, and ten). Of the remaining cases where the deterministic M/W model 
results were outside the boundary for a portion of the battle shown in the Total Force 
Level Trajectory plot, it was noted that the individual unit force level trajectories were 
within the established one standard deviation boundary in all but one case (case nine). 

Secondly, the results from the experiment show that both models force level 
trajectories displayed the expected mode of behavior for the number and type of 
equilibrium points which were present in all cases. This shows strong support for the 
research done by Professors Moose and Wozencraft. 

This experiments results are certainly not proof that the deterministic model is 
always a good approximation for the mean result of a stochastic model. However, the 


results do indicate they often may be a good approximation. 


DZ 


V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The results of this research showed that the M/W trajectory was a good 
approximation for the C/S trajectory in nine of ten cases. This is an indication that 
deterministic attrition models may be good approximations for the mean of stochastic 
attrition models. It has also shown that the modes of behavior for both models 
trajectories are very similar. The combination of these two results provides strong 
support that the M/W model is a good approximation for the C/S model. These 
results provide important support for the argument that a deterministic attrition model 
can be used to model the combat attrition process which is generally agreed upon as a 
stochastic process. 

This simulation provided a test of the Moose/Wozencraft model. The M/W 
model was tested utilizing four different force composition scenarios, for each of the 
four types of trajectory behavior modes. The results showed the M/W models results 
were good approximations for the individual unit level and the aggregate force levels of 
the C/S model. This Generalized Lanchester Equation type model is simple to use and 
understand. I[t allows the modeler, or military planner, to model an aggregate force of 
as many different type forces as required to model the real-world battlefield. The M/W 


model has real potential for future use by military analysts. 


B. COMMAND AND CONTROL APPLICATIONS 

A key decision in the Command and Control decision process is what force 
Structure is required to engage a proposed enemy. The problem is to determine the 
optimum force size and composition required to engage an enemy and win. Optimizing 
one’s force is a very complex problem, involving the right mixture of force types and 
the right number of each of the forces. An example would be when a commander may 
need a tank heavy unit for a specific battle in order to ensure that his force will win. 

A model such as the M/W model can serve as a very useful decision aid for the 
military planners and combat commanders to assist with optimization of forces on the 
battlefield. Given an accepted set of attrition-rate coefficients, the model can provide 
the capability of testing for the optimum mix and force size required to fight and win a 


battle. The use of a model of this type by maneuver units would require the 
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commanders and their staffs to ask the intelligence community for specific information 
on the enemies current force structure and resupply capability. This information could 
be used with the model to assist planners and commanders to decide whether to 
engage, and the optimum force compostion if the decision is made to do so. If not, 
they may chose to with-draw and fight on better terms. Models of this type would be 
beneficial to military planners responsible for developing OPLANS and CONPLANS 
also. It would allow them to hypothesize the expected enemy force structures and test 


proposed force structures to engage the enemy. 


C. RECOMMENDATIONS 

After completing this project it is apparent that a great deal of future research 
can be done with these two models to provide a more complete analysis. The results 
from this initial survey have identified areas for further testing which would provide a 
more complete understanding of the dynamics of the two models. The areas listed 
below would provide good topics for future work with these models: 

1)  Utihze the models to locate the equilibrium points and the separatrix for cases 
Where an unstable equilibrium point exits in the state space. Once they have 
been located, observe the two models results in the vicinity of the equilibrium 
point for several different inital force levels. The separatrix is a very 
interesting factor which should be observed closely. 

2) A detailed analysis of case nine of this experiment is recommended. Case nine 
was the one case which was run that did not behave similarly to the rest of the 
cases. At this time it is not Known why, and therefore deserves further 
research. 

3) One of the observations made from this experiment involved a difference in 
the attrition-rate of a force over time for the two models. This could be a 
result of the normalization factor used for the aimed fire attrition in the M/W 
model. This factor was used to ensure that an individual force was not 
allowed to fall below a zero level. This factor may cause the two models to 
have slightly different attrition-rates over time. This can be tested by 
removing the factor and changing the models computer program to stop if one 
force falls to a zero level. 


4) Test the model(s) against approved combat simulations. 
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Further research in this area may lead to an attrition model which is accepted as 
a useful model for modern combat. This could provide a useful Command and Control 
decision aid for use by commanders and their staffs to optimize force structures on the 
battlefield. 
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APPENDIX A 
EXEC FILE FOR C/S PROGRAM 


VFLODLIB 


INTAW DATA A 
TAW FTOQ6FOO] A 
TAW FTO7FOO1 A 
TATY DATA A 
X1Y1 DATA A 
X2Y2 DATA A 
X3Y3 DATA A 
STEP DATA A 
STATTF DATA 
STATF1l DATA 
STATFZ2 DATA 
STATF3 DATA 


PP py Ye 
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MY ANNANAAIAAAAANAINAAAINAAANIAIAAAANAIAAANAAAAAAAAARAAANRANARNNA 


WRITTEN 


Nige AG 
NXS ,YS(I 
NX,Y(I) 


AAA( I,J 
BBB(I,J 
COOLS 
Pep I,J 


APPENDIX B 
C/S COMPUTER PROGRAM 


PROGRAM LANCHESTER 
6/18/86 
TO TEST THE MOOSE/WOZENCRAFT THEORY 


NUMBER OF TYPES OF COMBATANTS <INPUT> 
) INITIAL NUMBER OF EACH TYPE OF COMBATANT 
NUMBER OF EACH TYPE OF COMBATANT ALIVE 


ATTRITION: X(I) LOSSES TO AREA FIRE 
ATTRITION: Y(J) SHOOTS X 
ALIRITION: Y(t) LOSSES TO AREA FIRE 
PU NiirenN:s ACI) SHOOTS *y 


NUMBER COMBATANT oe RESUPPLIED/TIMESTEP 
NUMBER COMBATANT Y(I) RESUPPLIED/TIMESTEP 


) STORED POSITION OF EACH COMBATANT 
K) STORED STATUS (LIVE OR DEAD) OF EACH COMBATANT 
STORED TYPE OF EACH COMBATANT 


J) PK: X SHOOTS Y TYPE (I,J) INTERACTION <INPUT> 
) AREA FIRE LOSS COEFFICIENTS FOR & <INPUT> 


LISTS OF RANDOM VARIANTS, CONTROLS WHO SHOOTS 
LISTS OF RANDOM VARIANTS, CONTROLS HITS 

LISTS OF RANDOM VARIANTS, CONTROLS SELF-LOSSES 
LISTS OF RANDOM VARIANTS, CONTROLS TARGET CHOICE 


PK: X(I) KILLS Y(J) THIS TIMESTEP 

NR OF . HIT BY oT, THIS TIMESTEP, AREA FIRE 
Pomomey o) Hele BY .A(l) THIS TIMESTEP, @AriED FIRE 
PROBABILITY OF SELF<-LOSS 


fou) LOSens [O ALIMEDGRIRES THIS, TIMESTEP 
K,Y(I) LOSSES TO AREA FIRE THIS TIMESTEP 
Meo Seer LOSsSskS THismiaiirSlEP 


MAXIMUM NUMBER OF TIMESTEPS <INPUT> 
POSSIBLE POSITIONS OF COMBATANTS <INPUT> 








IMPLICIT REAL*8 (D) 


CHARACTER*40 IHEAD 
CHARACTER*1 IANS,IYES 


DIMENSION 
DIMENSION 


DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 


Milena es, ) e747 (3222) NOX (3228) Bs) 
YPOS (3222) , YSTAT (3222) ,NOY(3222) ,NY(3 
aes} thee qCeee) 3 

YLOC(3), PKSY(3),NSSS(3 

RBE(3, 2}, XBR 3/3) 

PKX(3,3),PKY(3,3 

A 3/3), RUX(3,3) 0K (3.3) 
PKAK(3,3) , PKAY(3,3),KAX(3,3) ,KAY(3,3 
NDX(3) ,NDY(3) ,NXS(3) ,NYS(3) 

NDHX (3) ,NDHY(3 

NDAX(3) NDAY(3 

NDSX(3) NDSY(3 
AAA(3,3),AAAAV(3,3) , SAAAAV(3,3) 


a 


DIMENSION BBB(3,3),BBBAV(3,3),SBBBAV(3,3 
DIMENSION CCC(3,3),CCCAV(3,3).SCCCAV(3,3 
DIMENSION DDD(3,3),DDDAV(3,3) ,SDDDAV(3,3 

C 
DIMENSION RX(4222) RY (4222) , RN( 9688) 
DIMENSION Rx (4222), RY 
DIMENSION RXT(4222 RYE rope 

C 
DIMENSION SNTX (212), SNK1 (212) , SNH2 (212) , SNK3(212) 
DIMENSION SNTY(212),SNY1(212) ,SNY2(212) ,SNY3(212 
DIMENSION SSNIR( 212), SSNX1 (212), SsNKZ(212), SaNe3{ 212) 
DIMENSION SSNTY(212) .SSNY1(212) .SSNY2(212) .SSNY3(212 
DIMENSION ANTK (212), aNNL (212), ANK2(212), aNx stg tay 
DIMENSION ANTY(212) ,ANY1(212) ,ANY2(212) ,ANY3(212 
DIMENSION SDNTX(212) ,SDNX1(212) , SDNX2(212) , SDNX3(212 
DIMENSION SDNTY(212) ,SDNY1(212) ,SDNY2(212) , SDNY3(212 
DIMENSION PSNTX(212) , PSNX1(212) | PSNX2(212) , PSNX3(212 
DIMENSION PSNTY(212), PSNY1(212) , PSNY2(212) , PSNY3(212 
DIMENSION TSNTX(212). TSNX1(212) , TSNX2(212) , TSNX3(212 
DIMENSION TSNTY(212), TSNY1(212) , TSNY2(212) , TSNY3(212 
DIMENSION CSTEP(1003 

C 
DATA IYES/'Y!/ 

C 

C 


CALL EXCMS('CLRSCRN' ) 
READ (5,9400) IHEAD 
WRITE(4,9410 

9410 FORMAT(' DO YOU WANT TO CHANGE THIS HEADING? (Y/N) (A40)') 
WRITE(4, ae IHEAD 
READ (4,9405) IANS 

9405 FORMAT(A1) 
IF(IANS. MG .IYES) READ (4,9400) IHEAD 

9400 FORMAT (A40 

C 


e ASK USER FOR THE NUMBER OF REPLICATIONS DESIRED 
e 


WRITE 4 , 9490) 
9490 READ( ‘aia ae REPLICAS DO DESIRE? (MAX 100)') 
9495 FORMAT 14) 

WRITE(4,9480) 
9480 FORMAT (! DO YOU WANT FILES 6 AND 7 PRINTED? ') 

READ (4,9401) IANS 
9401 FORMAT Al) 

WRITE(6,9400) IHEAD 

WRITE 7 ,9400 IHEAD 

WRITE(8,9400) IHEAD 

WRITE(9,9400) IHEAD 

WRITE(10,9400) IHEAD 

WRITE(11,9400) IHEAD 

WRITE(12,9400) IHEAD 

WRITE(13,9400) IHEAD 

WRITE(14,9400) IHEAD 

WRITE(15,9400) IHEAD 

WRITE(16,9400) IHEAD 

READS, ea} NSTOP 

READ (5,9500 Nee NX,XLOC,PKSX,NRRR 

IF (IANS. E@.1YES)TH 

et ae 02) 

Ce 6,9501) NTYPX,NX,XLOC,PKSX,NRRR 


NDIF 

READ (5,9505) ( (PRHK(T, J) ,J=1,3),1=1,3) 
IF (IANS.EQ.IYES)TH 

WRITE (8, 9 ai) 

WRITE 6.9520) ((PKHX(I,J) ,J=1,3),I=1,3) 
NDI 

READ (5,9505) ( (KAE(T, J), J=1,3),151,3) 
IF (IANS .EQ.IYES )THE 
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QAAaAN 


AAN 


9510 
gorL 


Geren caaee 
WRITE(6,9520) ((XAF(I,J),J=1,3),I=1,3) 
ENDIF 


READ (5,9500) NTYPY,NY,YLOC,PKSY,NSSS 
IF(IANS.E . IYES ) THEN 

WRITE (6, 2 03 

WRITE(6,9501) NTYPY,NY,YLOC,PKSY,NSSS 


ENDIF 

READ (5,9505) EEE Dy DRv=1,3),1=1,3) 

IF(IANS.EQ.IYES) THEN 

Wee ree. 9 22) 

WRITE 6,9520) ((PKHY(I,J),J=1,3),I=1,3) 

ENDIF 

READ (5,9505) ((YAE(T, Peale ian 3) 

LF (TANS. E ,1YES)TH 

Wet TEye. 3 14 

WRITE 6,9520) ((YAF(I,J),J=1,3),I=1,3) 
F 


END 

FORMAT ee 931: } 

FORMAT(415 ,3F5.1,3F7.3,315 

FORMAT( | X TYPES NUMBER EACH LOCATIONS SELF-LOSS 

X RE-SUPPLY 

FORMAT( ! # Y TYPES NUMBER EACH LOCATIONS SELF-LOSS 

X RE-SUPPLY ') 

FORMAT aa 5) 

FORMAT(9F8.5,/) 

FORMAT(/,' KILL PROBABILITIES FOR FORCE X(I) SHOOTS BI) df 

X - 7! oe z are - PK(1,3 PK(2, 1) PK(2,2) PK(2,3)', 

X ' PK(3,1) PK(3,2) PK(3,3 

FORMAT (/,' KILL PROBABILITIES FOR cen Y(I) SHOOTS aE 

X eR aa ei Age 3 PK(2, 1) PR(D. "3 PK(2,3)', 
' PK(3,1) PK(3,2) PK(3,3 

FORMAT(/,' AREA FIRE COEFFICIENTS nace Y(I) SHOOTS x(J) ee 

ARCS 1) a )a'5) AA(1,3 (AA(Z2, 1) aa(2. 2m e,3)', 

! 3,1) AA(3,2) AA(3,3 

ECEE RVI AREA FIRE COEFFICIENTS als E X(I) SH S XC) Sf /, 

ee0l.1) €c(1,/2 (23)! , 


x | eg 't} aoa a} 4 ae cc (2, 1) cc(2, 2) ec 


INITIALIZE THE RANDOM NUMBER GENERATORS 


READ(5,9510) DX,DY,DN 
DSEEDX=DX 
DSEEDY=DY 
DSEEDN=DN 
IF(IANS. err . IYES ) THEN 
aeehe’ 4 
WRITE 6,9510) DX,DY,DN, DSEEDX, DSEEDY , DSEEDN 
Ee 8.3111) 
FORMAT ' - DOUBLE PRECISION RANDOM VARIANT SEEDS Ly 


INTEGER SEED VALUES') 
INITILAIZE THE STATISTICS VARIABLES 
DO 10 I=1, 212 

SNTX(TI )=0 
SNX1 (I )=0 
SNX2(I )=0 
SNX3 (1 )=0 
SNTY(I)=0 
SNY1 (I )=0 
SNY2(1I)=0 
SNY3(I}=0 
SSNTX =0 
SSNX1 
SSNX2 
SSNX3 
SSNTY 


ob 


NANAK 


10 


7 


PSNK1 
PSNX2 
PSNX3 
PSNTY 
PSNY1 
PSNY2 
PSNY3 
TSNTX 
TSNA1L 
TSNX2 
TSNX3 
TSNTY 
TSNY1L 
TSNY2 
TSNY3 


CONTINUE 
DO 15 I=1,3 
NYS tans 

X(I 


NXS (I 


15 CONTINUE 
MINST=99 
AVSTEP=0.0 
SDSTEP=0.0 
SSTEP=0 .0 
SSSTEP=0.0 


DO 17 I=1,100 
CSTEP(I)=0.0 


CONTINUE 
BEGIN REPLICATIONS 


DO 4000 N=1,L 


NSTEP=0 
NTX=0 
NTY=0 
20 I=1,3 
I)=NYS(1I 
I)=NXS (I 


DO 
NY 
NA 
IF 


IF(I.LE.NTYPY) 


I.LE.NTYPX) 





E 


OOODWDVOWVODOVOMOOVOODO0O9O0O00000090 


=N 


THEN 


NTX=NTX+NX (TI) 
ENDIF 


THEN 


NTY=NTY+NY (TI) 
ENDIF 

20 CONTINUE 
NTXS=NTX 
NTYS=NTY 
NXSTOP=50 
NYSTOP=50 


60 


QAAaAN 


AAaAN 


AAAAANANAAN 


INITIALIZE THE TWO FORCES 


IMAX=0 
IMIN=1 
IP=0 


DO 200 J=1,NTYPX 
IMAX= IMAX+NX(J) 
DO oe I=IMIN, IMAX 


[Pai et 


SPOS(IP)= KLOC (J) 


XSTAT(IP)=1 

NOX(IP)=J 
100 CONTINUE 

IMIN=IMAX+1 
200 CONTINUE 


IMAX=0 
IMIN=1 
IP=0 


DO 400 J=1,NTYPY 
IMAX=IMAX+NY (J) 
DO 300 I=IMIN, IMAX 


IP=IP+1 


YPOS(IP) Foes 


YSTAT(IP)=1. 

NOY(IP)=J 
300 CONTINUE - 

IMIN=IMAX+1 
400 CONTINUE 

DO 500 I=1,3 

NDX(I)= 


Co 

e8) 

Co 

1 
< Oo 
© 

i 

Cy 
iu ius s 


DDDAV(L,J 
700 CONTINUE 


SETUP 


OOOOWUW 
OOOO 


A FILE FOR FUTURE PLOTTING 


ARITE(?. 79900) . LYES ) THEN 


WRITE 
9700 FORMAT 


fe RtKS NTX NX1l NX2 NX3 DX1 DX2 DX3', 


Neco,  NYL NyZ eels DYL DYZ DY3') 


x 
WRITE 7 9710) NSTEP,NTXS ,NTX,NX,NDX,NTYS ,NTY,NY,NDY 


9710 FORMAT 
oe | aaiiian | # 


H(A 
WRITE(6,97115 
‘NTXS NTX NX1 NX2 NX3 RX1 RK2 RX3!, 


2(515,314),/) 


NiiooeNy Weetine = NYS RY) RY2Z RY3/) 


NTYS 
WRITE (6, 9710) NSTEP ,NTAS ,NTX,NA,NRRR,NTYS,NTY,NY,NSSS 
N 


NTF |. 
NSX,Y¥ 

I a 
NTARX,Y 
NCX,Y 
Je 


TOTAL NUMBER OF LIVE COMBATANTS 
NUMBER OF COMBATANT WHO IS FIRING 
TYPE OF COMBATANT FIRING 

COUNT OF TARGETS ATTACKED 

NUMBER OF TARGET UNDER FIRE 

TYPE OF TARGET UNDER FIRE 


THIS IS THE BEGINNING OF THE COMBAT CYCLE 


1000 CONTINUE 


NSTEP=NSTEP+1 


NTT=NTX+NTY 


NGT=1.5*REAL (NTT) 
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AAD 


AANA 


AARAAAAAAANAANARANAARAANANNN 


Na) Seeeeeee 
NGY=1.3*REAL(NTY 


MAKE THE LISTS OF RANDOM NUMBERS 


CALL GGUBS Bee tae 
CALL GGUBS (DSEEDX ,NGA, RA 

CALL GGUBS (DSEEDA,NG&, ce 
CALL GGUBS(DSEEDX ,NGA&, RAT 
CALL GGUBS(DSEEDY,NGY,RY) 

CALL GGUBS (DSEEDY ,NGY, fet 
CALL GGUBS ( DSEEDY , NGY , RYT 


ZERO STORAGE FOR THIS TIMESTEP 


NSX=0 
NTARX=0 
NSY=0 


NDAY (TI )=0 
NDSX(T )=0 
NDSY (I )=0 
DO 1020 J=1,3 
KAX(I,J)=0 
KAY (I,J )=0 
KHX(I,J )=0 
KHY(I,J)=0 
1020 CONTINUE 
1040 CONTINUE 
DO 1060 I 
DO 1060 
PKAX(I, $4 
PKAY (I,J 
1060 CONTINUE 
1080 CONTINUE 
DO 1090 J=1,3 
DO 1090 I=1,3 
ST ety Gace: ae 
PKY(I,J)=PKAY(1I,J)+PKHY(1I,J 
1090 CONTINUE 


=1,3 
J=1,3 
mre ie 
=YAF(I,J)*NTX 


DO AIMED AND AREA FIRE INTERACTIONS FOR THIS TIMESTEP 
CHOOSE THE COMBATANT WHO SHOOTS, USING RN(I) 
MAKE COMBATANT NUMBER (COUNT), NSX,Y 


COMBATANT MUST BE ALIVE 


DETERMINE COMBATANT TYPE, LX,Y 
COMBATANT MUST NOT BE A SELF-LOSS 


INCREASE THE TARGET COUNT, NTARX,Y 
MAKE NUMBER OF TARGET, NCX,Y, WITH RA&AT,YT 


DETERMINE TARGET TYPE, L&,Y 


COMPARE : Egg ert > RX,Y, TARGET IS HIT 
jis < RX,Y, MISS 
PKHX,Y(LX,LY) > RX,Y, HIT BY ATMED ae 
< RX,Y, HIT BY AIMED FIRE 


Ay tN UA it 


PRA 2 Ce 
1100 CONTINUE 
DO 1180 M=1,NGT 
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AMO 


ENDI 
SE 


IF(RN(M).LE.(REAL(NTX)/NTT)) THEN 


X SHOOTS 
NSX=NSX+t1 
i RY GO TO 1120 
IF (XSTAT(NSX).EQ.0.0) GO TO 1180 


LX=NQX (NSX 

IF(PKSX(LX).LT.RXS(NSX)) THEN 
NTARY=NTARY+1 
NCY=RYT (NTARY) *NTY+1 
LY=N et 
IF(PKX(LX,LY).GE.RX(NSX)) THEN 


KILL TARGET NUMBER NCY 


IF(YSTAT(NCY).GT.0.0) THEN 
YSTAT(NCY )=0.0 
IF (PKHX(LX,LY).GE.RX(NSX)) THEN 
NDHY (LY )=NDHY (LY) +1 
pe LE BREA (LA LY) +2 


NDAY (LY )=NDAY(LY)+1 
KAX (LX, LY)=KAX(LX,LY)+1 
ENDIF 
ENDIF 


ENDIF 
GO TO 1180 
XSTAT(NSX)=0.0 


NDSX(LX)=NDSX(LX)+1 
GO TO 1180 
F 


ONO LS.- 


NSY=NSY+1 

ees? -S1 NTY) GO TO 1150 

IF(YSTAT(NSY) .EQ.0.0) GO TO 1180 

LYONOY (NSY 

IF(PRSY(LY).LT.RYS(NSY)) THEN 
NTARX=NTARK+1 
NCK=RXT (NTARX) ANTX+1 
a) 
IF(PRY(LY,LX).GE.RY(NSY)) THEN 


KILL TARGET NUMBER NCX 


IF (XSTAT(NCX) .GT.0.0) THEN 
XSTAT(NCX)=0.0 
IF (PKHY(LY,LX).GE.RY(NSY)) THEN 
NDHX(LX)=NDHX(LX)+1 
KHY (LY, LX)=KHY (LY, LX) +1 


NDAX(LX)=NDAX(LX)+1 
BO) HRCE Le)? 


ENDIF 
ENDIF 
GO TO 1180 
SE 
YSTAT(NSY)=0.0 


NDSY(LY)=NDSY(LY)+1 
GO TO 1180 


ALL TYPE X COMBATANTS HAVE FINISHED FIRING 
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1120 CONTINUE 
KSTOPX=1 
GO TORE 


ALL TYPE Y COMBATANTS HAVE FINISHED FIRING 


1150 CONTINUE 
KSTOPY=1 


EVERY ONE HAS FINISHED FIRING 


1170 CONTINUE 
IF(KSTOPX.EQ.1 .AND. KSTOPY.EQ.1) GO TO 1200 
1180 CONTINUE 


CALCULATE THE KILL RATIOS 


1200 CONTINUE 
DOVEZIO J=1, Nex 
DORI Z 50. L=1 NIP Y 
NUMD= NX(J) *NY(I) 
IF (NUMD.GE.1) THEN 
DDD(J, I )=REAL(KHX(J,1)*NTY) /REAL(NUMD ) 


E 
DDD(J,1)=0.0 
ENDIF 
DDDAV(J, iS eee i) DOP) 
NUMB=NY (I) *NX(J le 
IF (NUMB.GE.1) T 
BBB(I,J)= =KEAL (KHY (I, J)*NTX) /REAL(NUMB ) 


BBB(I,J)=0.0 
ENDIF 
BBBAV(I, a3 [BBBAV(T, J)+BBB(I,J) 
NUMC=NX(J)*NY(I 
IF (NUNC. GE.1) THEN 
CCC(J,1)=REAL(KAX(J,1))/REAL(NUMC) 


E 
ccc(J,1)=0.0 


ENDIF 
CCCAV( J, 1s ec ean i) + CCG) am 
NUMA= NY (I ANX(J) 


IF (NUMA.GE.1) THEN 
AAA(I, 7 =REAL(KAY(1I,J))/REAL (NUMA) 


SE 
AAA(I,J)=0.0 
DIF 


AAAAV(I,J)=AAAAV(I,J)+AAA(I,J) 
1250 CONTINUE 


TOTAL KILLED BY ALL MECHANISMS 
ee 
NDY( J )=NDHY(J)+NDAY(J)4+NDSY(J 
1290 CONTINUE 


IF(IANS.EQ.IYES)TH 
WRITE (6, 9660) NSTEP. ((KHX(I,J) ,J=1,3),1=1,3) 


ENDIF 
9660 FORMAT(/, TIMESTEP', 14,/, 
x AIMED mee ‘KILLS OF Y"S B 


Y & FORGE 
X MT EF ) a 3} Dba 33 D( 2 pip eee ee ee (2,3) '% 
x D on Bs 2 D(3,3)) 7/7 
x 
IF(IANS. Taee I¥ES) 
WEEE, 9665) CCKHY (L, J) ,J=1,3),1I=1,3) 


ENDIF 
9665 FORMAT(10x, | AIMED FIRE KILLS OF X"S BY Y a 
X //," Ba1,1) B(1,2) B(1,3) B(2, ays *3 (2,2) B(2,3)', 
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QaAaAN 
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9662 


9668 


B601 


Boo0 


1500 


1510 


9670 


160 


1620 


0 


X tienen ( 3,2) B(3,3)',/,1%, 
918, 


x 
IF( IANS. nae - LYES) 
WRITE 6, 9662) OURAN CL, J), J=1,3),1I=1,3) 


END 
FORMAT (10X | AREA FIRE KILLS OF Y"S BY X FORCE 7 
x /},' (1,1 c(t, 2} cf, 2 et) "elae)  c(2,3)', 

1 Sih 3,2 3°3)' /, 1x, 


(IANS. “50. .IYES)TH 
WRITE 6, 9668) i ere. J),J=1,3),1=1,3) 
FORMAT(10X,' AREA FIRE KILLS OF X"S BY Y FORCE 7 
X pipes Bit} ACL, 2) ve ao a(2 t) Comma 2.3)! 
1 3,1 3,2) A(3,3)',/ 


IF (IANS .EQ.IYES)THEN 


, 


vl Sy 


WRITE(6,9700 

WRITE(6,9710) NSTEP,NTXS ,NTX,NX,NDX,NTYS ,NTY,NY,NDY 
WRITE(7,9651 

ENDIF 

FORMAT(' D,B,C,A ves ROWS IN DESCENDING ORDER ') 
Pen 29:1 .IYES)TH EN 

WRITE NS Raat sd) p»g=1l,3),L=155 
WRITE 3650 NSTEP, hal) :J=1,3),.=1.,5 
WRITE +9680 NSTEP, ((KAK(I,J),J=1,3),1=1,3 
WRITE(7 ,9650 NSTEP, KAY(I,J),J=1,3),1=1,3 
ENDIF . 

FORMAT (14,918) 


CALCULATE TOTAL LOSSES FOR ALL INTERACTIONS 


CONTINUE 
NDTX=0 
NDTY=0 
DO 1510 J=1,3 
NDTY<NDTYNDY(3} 
NDTY=NDTY+NDY (J 
CONTINUE 
IF(IANS. Ee . LYES ) THEN 
wEtTE Ce, 9670) NDTX,NRRR,NDHX,NDAX,NDSX, 
x NDTY, ‘NSSS, ‘NDHY, ‘NDAY, ‘NDSY 
WRITE(7,9670) NDTX, NRRR, NDHX, NDAX, NDSX, 
N \NSSS, ‘NDHY, ‘NDAY, NDSY 
FORMAT(' TOTAL X FORCE LOSS WAS |, TAY RESUPPLY WAS:',314,/, 


X 1 AIMED FIRE LOSSES:' , 314, 
x 1 AREA FIRE LOSSES:!' ,314,/, 
X 1 SELF LOSSES:' , 314, /, 
X ' TOTAL Y FORCE LOSS WAS',14,', RESUPPLY WAS:',314,/, 
X 1 AIMED FIRE LOSSES:!' ,314,/. 
X 1 AREA FIRE LOSSES:!',314./, 
X 1 SELF LOSSES:! 314) 

ENDIF 

ADD IN RESUPPLIED UNITS 
CONTINUE 
DO 1650 K=1, 


IF(K.LE. NTYDX »~AND. NRRR(K).NE.0) THEN 
NADDX=NRRR (K 
NX(K)=NX(K)+NADDX 
DO 1620 I=1,NADDX 
NTX=NTA+1 
Meare =XLOC(K) 


N 
RS TAT NTXOE =1.0 
CONTINUE 


ENDIF 
IF(K.LE.NTYPY .AND. NSSS(K).NE.0O) THEN 
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NADDY=NSSS(K) 

NY (K)=NY(K)+NADDY 
DO 1640 I=1,NADDY 
NTY=NTY+1 
2 =YLOC(K) 


N 
vSTAT (NYY): =1.0 
1640 CONTINUE 
ENDIF 
1650 CONTINUE 


CLOSE OUT DESTROYED UNITS 


1700 CONTINUE 
DOma LO Jaks 
NA( J )=0 
NY (J )=0 
1710 CONTINUE 
TX=0 


DO 1740 I=1,NTX 
TE ee .NE.O. 0) THEN 


ie Lt 

NOX(TX)= NBs) 
XPOS (IX)=XPOS(1I) 
XSTAT(IX)= RSTAT(I) 


K=NOX(I 
NX(R)=NX(K)+1 
ENDIF 
1740 CONTINUE 
NTX=1X 
IY=0 
DO 1780 I=1,NTY 
IF(YSTAT(1I).NE.O. ° THEN 
he i¥+] 
NOY (TY) NOY (I 
YPOS (IY)=Y¥POS( (1) 
YSTAT - =YSTAT(I) 


Nv(R) =NY(K)+1 
Dir 


1780 CONTINUE 
NTY=TIY 


PRINT STATUS FOR EACH TIME STEP IF REQUESTED 


IF (IANS .EQ.IYES) THEN 
ee bade oe NSTEP ,NTAS ,NTX,N&X,NDK,NTYS ,NTY ,NY,NDY 


WRITE(6,9710) NSTEP,NTXS ,NTX,NX,NRRR,NTYS,NTY,NY,NSSS 
END IF 
ACCUMILATE STATISTICAL VARIABLES FOR EACH TIME STEP 
I=NSTEP 
ONTX (I )=SNTA(I)+NTX 
SONAL (I )=SNAL(TI)+NX(1 
SNX2(T )=SNX2(1)+NX(Z 
SNX3 (I )=SNK3(1)+NX(3 
SNTY (I )=SNTY (I )+NTY 
SNY1(I)=SNY1(1I)+NY(1 
SNY2(I)=SNY2(1I)+NY(2 
SNY3 (I )=SNY3(1)+NY(3 
SSNTX(1I)=SSNTX(1I)+NTX**2 
SSNX1(1I)=SSNX1 (1 )+NX(1)**2 
SSNX2 (I )=SSNK2(1I)+NX(2)**2 
SSNX3 (I )=SSNX3 (I )+NX(3 ) **2 
SSNTY(I)=SSNTY(I )+NTY**2 
SSNY1(I)=SSNY1(1)+NY(1)**2 
SSNY2 (I )=SSNY2(1)+NY(2) **2 
SSNY3(1)=SSNY3(1)+NY(3)**2 
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TERMINATE IF FORCES ARE ATTRITTED BELOW SPECIFIED LEVEL 
IF(NTX.LE.NXSTOP .OR. NTY.LE.NYSTOP) GO TO 9120 
TERMINATE ON NUMBER OF CYCLES 


IF(NSTEP.LT.NSTOP) GO TO 1000 
9120 CONTINUE 
STEPS=REAL (NSTEP) 


C 
C CALCULATE THE AVERAGE ATTRITION COEFFICIENTS PER REPLICATION 
ce 


DO 9100 J=1,3 
DO 9100 I=1,3 
AAAAV(I,J)=AAAAV(I,J)/STEPS 
BBBAV(I,J)=BBBAV r’ .J)/STEPS 
@GCAy( 1, Jy=eecay L, J /STEPS 
DDDAV(I,J)=DDDAV(I,J)/STEPS 
SAAAAV(I,J SS AAAAV( I,J)+AAAAV(I, 
SBBBAV(I,J)=SBBBAV(I ,J)+BBBAV I, 
Le 


AAD ANN 


SCCCAV(T,J)=SCCCAV(T,J)+CCCAV 

SDDDAV(I,J)=SDDDAV(1I,J)+DDDAV 
9100 CONTINUE 

CSTEP(N)=NSTEP 


C 
C ACCUMILATE THE STATISTICAL VARIABLES FOR THE NUMBER OF TIME STEPS 
C 
SSTEP= SSTEP+CSTEP(N) 
: SSSTEP=SSSTEP+(CSTEP(N)**2) 
: WRITE THE NUMBER OF REPLICATIONS AND THE NUMBER OF STEPS TO FILE 


sts Bai iN ,CSTEP(N) 
8510 FORMAT(I3,F4 
ie 


e KEEP TRACK OF MINIMUM NUMBER OF STEPS FOR ALL REPLICATIONS 
C 


IF(MINST.GT.CSTEP(N) ) THEN 
MINST=CSTEP(N) 


ENDIF 
4000 CONTINUE 
c 
CALCULATE AVERAGE/STANDARD DEV.OF FORCE LEVELS FOR EACH TIME STEP 
B=REAL(L) 
RD=1.0/B 
RDD=1.0/(B*(B-1) ) 


DO 4400 I=1, MINST 
ANTX(I)=SNTX(I)*RD 


ANX1(I)=SNX1(I)*RD 

ANX2(I)=SNX2(I)*RD 

ANX3(I)=SNX3(I)*RD 

ANTY(I)=SNTY(I)*RD 

ANY1(I)=SNY1(I)*RD 

ANY2(1T)=SNY2(1I)*RD 

ANY3(I )=SNY3(1I)*RD 

SDNT =SORT (ABS ( (B*SSNTX(I) -SNTX(I)**2)*RDD 
SDNX1(1I)=SORT (ABS ( (BXSSNX1 (I) -SNX1(1I)**2)*RDD 
SDNX2(I)=SORT (ABS ( (B*SSNX2(I)-SNX2(I)**2)*RDD 
SDNX3(I)=SORT(ABS ( (B*SSNX3 (I) -SNX3(I)**2)*RDD 
SDNTY(I)=SORT(ABS( (BXSSNTY(I )-SNTY(I)**2)*RDD 
SDNY1(1I)=SQORT(ABS( (BXSSNY1(1)-SNY1(1I)**2)*RDD 
SDNY2(I )=SORT(ABS( (BXSSNY2(I)-SNY2(I)**2)*RDD 
SDNY3(I)=SORT(ABS ( (B*SSNY3(1I)-SNY3(I)**2)*RDD) 


CALCULATE THE VALUES FOR PLUS AND MINUS ONE STANDARD DEVIATION 
FROM THE MEAN VALUES FOR PLOTTING 


AAAN 
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PSNTX(1I)=ANTX(1)+SDNTX(I 
PSNX1(1)=ANX1(1)+SDNXL(I 
PSNX2(1)=ANK2(1)+SDNX2(T 
PSNX3(1)=ANX3(1)+SDNX3(I 
PSNTY(1I)=ANTY(1)+SDNTY(I 
PSNY1(1I)=ANY1(1)+SDNY1(T 
PSNY2(1 )=ANY2(1)+SDNY2(T) 
PSNY3(1)=ANY3(I +SDNY3(T) 
TSNTX(1I)=ANTX(1)- (SDNTX(T) ) 
TSNX1(1I)=ANX1(1)-SDNX1(I 
TSNX2(1)=ANX2(1)-SDNK2(I 
TSNX3(1)=ANX3(1)-SDNX3(I 
TSNTY(1I)=ANTY (1)-SDNTY(I 
TSNY1(1)=ANY1/1I)-SDNY1(I 
TSNY2 (1) =ANE2 1) -SDNE2 I 
TSNY3(1)=ANY3(1)-SDNY3(T 


c 
CALCULATE THE AVERAGE ATTRITION COEFFICIENTS OVER ALL REPLICATIONS 


DO 9200 J=1,3 
DO 9200 K=1,3 
AAAAV(K, J )=SAAAAV(K, 
BBBAV(K,J)=SBBBAV(K, 
CCCAV(K, J )=SCCCAV(K, 
DDDAV(K,J)=SDDDAV(K, 
ee CONTINUE 


EC PRINT STATISTICS TO FILES 8,9,10,11 BORSEEotiaNG 


nee B80) I ,ANTX(1) ,ANZY(1),,PSNTX(1) , PSNDY( 1). TSNIxX( 
XTSN 
eo 9,8100) I,ANX1(1I),ANY1(1I),PSNX1(I),PSNY1(1I),TSNX1(I), 
SN i 
neta Ry I,ANX2(1I),ANY2(1I),PSNX2(1I),PSNY2(1I) ,TSNX2(TI), 
SN 
WRITE(11,8100) I,ANX3(1I),ANY3(1I),PSNX3(1I),PSNY3(1),TSNX3(I), 
onto 
8100 FORMAT(I3,6(F8.2,1X)) 


Cc 
; WRITE AVG'S AND STANDARD DEVIATIONS TO FILE FOR TABLES 


J) *RD 
J) *RD 
J)*RD 
J)*RD 


WRITE(13,8150)I,ANTX(I) ,SDNTX(I) ,ANTY(I) ,SDNTY(I 
WRITE(14,8150)1I,ANX1(1I),SDNX1(1I),ANY1(I),SDNY1(I 
WRITE(15,8150)1I,ANX2(1I),SDNX2(1I),ANY2(I) ,SDNY2(I 
WRITE(16,8150)1I,ANX3(1I),SDNX3(1I) ,ANY3(1I),SDNY3(I1 


8150 FORMAT(1I3,4(F8.2,1X)) 
4400 CONTINUE 


c 
C CALCULTE THE AVG AND STANDARD DEVIATION FOR THE # OF NSTEPS PER RUN 
C 


AVSTEP=SSTEP*RD 
SDSTEP=SQRT (ABS ( (B*SSSTEP-SSTEP**2) *RDD) ) 
WBE TS Soa cgy ho. -AVSTER, SBOREE 

8500 FORMAT(1I3,2F7.2) 


EXPLAIN CAUSE OF TERMINATION 


9000 CONTINUE 
WRITE(7,9 

9656 FORMAT('! 
WRITE 
WRITE 


5 

9655 FORMAT(' A 
WRITE 2 

Dia) 


WRIT 
9690 FORMAT(/,' TIMESTEP', 14,/, 
X 10X,' AVERAGE ATTRITION CO 
X | era's com GE 
X ' GCs, 1) Ceiter mee 


AAD 


Bye Ronse Uh 
eens = 
BBBAV(I,J),J= 
o) 

eee as 


KX EtOs) "ee 


FOR X FORCE 
1 eZ C(Z, Soa 


K SP, Ore .5,/) 
WRITE (7 ,9655 See5V([ J) ,J-1,3),1=1,3 
WRITE(7,9655 DDDAV(I,J),J=1,3),1=1,3 
WRITE(6,9698 CCGAVC lo J=-1,3),1=1,3 
WRITE(6,9699 DDDAV(I,J) ,J=1,3),1I=1,3 
9695 FORMAT(10xX, 
x ' AVERAGE ATTRITION COEFFICIENTS FOR 
x oe nate 3} ee AA(2,1) 
x wee AACS, 2) BAC3,3)',/,1%, 
SF Sion 


x 
poo] BO ee 


' AVERAGE ATTRITION COEFFICIENTS FOR 
X ema mara | AE BB(2,1) 
x a Bee) ebb. 3s,2) BB(3,3)' ,/ aLx, 
FS.5, 


x 
9698 FORMAT(10X, 
x ' AVERAGE ATTRITION COEFFICIENTS FOR 


x pea Bera’ s} Bela s COC Zs) 
x meer iemec (sa 2) CC(3,3)",/7,1%, 

SES 5, / 

9699 FORMAT(10X, 
x ' AVERAGE ATTRITION COEFFICIENTS FOR 
x eoratd Bess Spa's DOKZ, 1) 
x PDK Ss, )) wees, 2) DD(3;73)',/,1%, 
x SS 25) 
Pp 


IF(NTX.LT.NXSTOP) THEN 
WRITE $8 36853 NSTEP ,NTA,NTAS 
WRITE (7,9620) NSTEP,NTX,NTXS 
END 
IF(NTY.LT.NYSTOP) THEN 
WRITE Me oeaat NS CER Ne Na 
ete 7, 9630) SiStEP NITY ,NTYS 


IF (NSTEP .EQ. NSTOP) THEN 
WRITE ie ceaas NSTOP 
WRITE (7,9640) NSTOP 


FORCE 
(2,2) 


a 


AA(2,3)' 


Y FORCE 
Be(2, 2)eeb«2, 3)! 


XK FORCE 
CeCZna)eee( 2,3)" 


K FORCE 
BDZ pe D 2, 3)' 


ENDIF 
9620 FORMAT(' STOPPED IN STEP ',1I5,' BY FORCE X',/, 
x ' CURRENT X STRENGTH =',15,', INITIAL X STRENGTH! ,1I5,//) 
9630 FORMAT(' STOPPED IN STEP ',I5,' BY FORCE Y',/, 
x ' CURRENT Y STRENGTH =',15,', INITIAL Y STRENGTH! ,I5,//) 
9640 FORMAT(' STOPPED BY PROGRAM AT NSTOP = ',15,//) 


END 
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APPENDIX C 
EXEC FILE FOR M/W PROGRAM 


&TRACE ON 

GLOBAL TATLIB VFORTLIB NONIMSE IMSESF CMSETE 
GLOBAL LOADLIB VFLODLIB 
FILEDEE 4) 2ERM 

FILEDEF 5 DISK INMWM DATA A 
FILEDEP 6 DISK Miia sTE teat en 
FILEDEF 8 DISK TOTAL DATA A 
FILEDEF 9 DISK FOR1 DATA A 
FILEDEF 10 DISK FOR2 DATA A 
FILEDEF 11 DISK FOR3 DATA A 
LOAD MWM 

START MAIN 

CP MSG 0731P JOB MWM COMPLETED 
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APPENDIX D 


M/W COMPUTER PROGRAM 


PROGRAM WRITTEN TO RUN THE GENERALIZED LANCHESTER 
EQUATION MODEL DESIGNED BY 
PAUL H. MOOSE AND JACK M. WOZENCRAFT 


1/08/87 


WRITTEN TO TEST THE MOOSE/WOZENCRAFT THEORY 


NTYPX 
Nex PY 
mes’ 
oN J 
SINK, ¥Y 
iNA 5 2 


AAA(J,I 
BBB i I 
SEC L; a] 
DDD I,J 


SAAA 
SBBB 
ScCCe 
SDDD 





NUMBER OF TYPES OF XK COMBATANTS <INPUT> 
NUMBER OF TYPES OF Y COMBATANTS <INPUT> 
NUMBER OF EACH TYPE OF X COMBATANTS 





NUMBER OF EACH TYPE OF Y COMBATANTS 
VARIBLE TO SUM TOTAL FORCE LEVEL PER TIMESTEP 
TOTAL FORCE LEVEL PER TIMESTEP 
ATTRITION: a AREA FIRE LOSSES DUE TO Y(J) 
ATTRITION: X AIMEDSEITRE LOSSES DUE TO Y(J 
ATTRITION: Y 5 ARERA FIRE LOSSES DUE TO X(I 
ATTRITION: Y(J) AIMED FIRE LOSSES DUE TO X(I 
ATTRITION: A(I) SUM OF THE AREA LOSSES 
ATTRITION: X(I) SUM OF THE AIMED FIRE LOSSES 
ATTRITION: Y(J) SUM OF THE AREA FIRE LOSSES 
ATTRITION: Y(J) SUM OF THE AIMED FIRE LOSSES 
NUMBER COMBATANT Ca RESUPPLIED/TIMESTEP 
NUMBER COMBATANT Y(J) RESUPPLIED/TIMESTEP 
SELF~ATTRITION COEFFICIENT FOR X FORCE 
SELF-ATTRITION COEFFICIENT FOR Y FORCE 

KCl eLOssihS TO AIMED FIRE THIS TIMESTEP 

Y¥(J) LOSSES TO AIMED FIRE THIS TIMESTEP 

BH LOSSES TO AREA FIRE THIS TIMESTEP 

Y(J) LOSSES TO AREA FIRE THIS TIMESTEP 

2(J) SELF LOSSES THIS TIMESTEP 

TOTAL CHANGE IN as FORCE LEVEL THIS TIME STEP 


7. SELF LOSSES THIS TIMESTEP 


TOTAL CHANGE IN Y(J) FORCE LEVEL THIS TIME STEP 


MAXIMUM NUMBER OF TIMESTEPS <INPUT> 
X FORCE LEVEL DEFINED AS THE BREAK PT. 
Y FORCE LEVEL DEFINED AS THE BREAK PT. 


IMPLICIT REAL*8 (D) 


CHARACTER*40 IHEAD 
CHARACTER*1 


DIMENSION TNX( 
DIMENSION NRRR 
DIMENSION AMLX 
DIMENSION ARLX& 


IANS , IYES 
3) 


TNY(3) 
»NSSS 
, AMLY 
, ARLY 


3 
) 
3 








3 
3 
3 


TI 


CY¥O, 0107 


eeke! 


DIMENSION SLX(3) 


SLY (3) 
DIMENSION TNDX(3), TNDY(3) 
DIMENSION AAA(3,3 v GAARA 3 
DIMENSION BBB(3,3),SBBB(3 
DIMENSION ccc(3,3),SCCC(3 
DIMENSION DDD(3,3),SDDD(3 
DIMENSION UUU(3$,VVV(3) 
DATA IYES/'Y'/ 
CALL EXCMS('CLRSCRN' ) 
READ (5,9400) IHEAD 
WRITE(4, 9410 


9410 FORMAT(' DO YOU WANT TO CHANGE THIS HEADING? (Y/N) (A40)') 


9405 


9400 
9401 


WRITE(4, Biss IHEAD 

READ (4,9405) IANS 
FORMAT (A1) 

IF LANS - ROMA) READ (4,9400) IHEAD 
FORMAT (A40 
FORMAT (Al ) 

WRITE(6,9400) IHEAD 

WRITE(8,9400) IHEAD 

WRITE(9,9400) IHEAD 

WRITE 10,3400) IHEAD 

WRITE(11,9400) IHEAD 

READ (5,9500) NSTOP 

READ (5,9500) NTYPX,TNX,UUU,NRRR 
WRITE(6 ,9502 

WRITE(6,9501) NTYPX,TNX,UUU,NRRR 

READ (5,9505) ((DDD(I,J},J=1,3),I=1,3) 
WRITE(6,9521 

WRITE(6 ,9520 ((DDD(T, J}, ie 3), /I=1, 3) 
READ (5,9505 CCC(T, 3) ,J=1,3)ei=rs 
WRITE (6 ,9523 

WRITE(6,9520) ((CCC(I,J),J=1,3),I=1,3) 
READ (5,9500) NTYPY, TNY,VVV,NSSS 
WRITE (6 ,9503 

WRITE(6,9501) NTYPY,TNY,VVV,NSSS 

READ (5,9505) ((BBB(1L,J),J=1,3) ,I=1,3) 
WRITE(6, 9522 

WRITE(6, 9520 { (BBB(T, J) Jot, 3)4T81,3} 
READ (5,9505 AAA( I Se O=1 cee S 
WRITE(6 ,9524 

WRITE(6,9520) ((AAA(I, J) Jad), pei 3) 
FORMAT(15,3F5.0,3F7.3,315 
FORMAT(1I5,3F7.2,3F7.3,315 
FORMAT ( 


FORMAT 3F8. 5) 

FORMAT (9F8.5,/) 

FORMAT (/," KILL PROBABILITIES. FOR FORCE X(I) SHOOTS 
x Bh’ 4 ra Mt: ie 2 DD(2, 1) DD( 252 
x ' DD(3.1) DDi2,2) DDS 

FORMAT(/,' KILL PROBABILITIES SOR FOR OOTS 
x BEI) ah 2) co a BB(2, iy) BBs. 2) 
X ' BB(3,1) BB(3,2) BB(3,3 

FORMAT(/,' AREA FIRE COEFFICIENTS FORCE Y(I) SHOOTS 
X VAR(L 1) AA 1,3) AA au AA(2 1) sacle 
X AA(3,1) AA(3,2) AA(3,3)') 

FORMAT(/,"- AREA FIRE COEFFICIENTS FORCE X(I) SHOOTS 
x 1 nae 1} ae c(h Cle. 1) Cala?) 
iz ' ¢c(3,1) €e(3,2) CC(s ape 


INITIALIZE THE INTIAL FORCES LEVELS 


# X TYPES NUMBER BACH SELF=EO55 
x RE=SUPPLY 

FORMAT( ' iF Yor weES NUMBER EACH SSE LE Snes 5 
x RE=SUrRL " 


We 


? 


? 


INTA=0 
INTY=0 
DO=rOOmrat. 3 
tt 
INTY=INTY+INY (I 
1100 CONTINUE 
TINA=INTA 
TINY=INTY 
XSTOP=50.0 
YST@E—50.0 
MINST=99 


SETUP A*PILE FOR TRACING EACH STEPS RESULTS 


Meee 

9711 FORMAT # TES TNX TNX2Z INX3 RX1 RK2 RX3', 

K SU tieeeriice tNYT TNY2 TNY3 RY1l RYZ RY3') 

2) NSTEP,INTX, TINX,TNX,NRRR,INTY,TINY,TNY,NSSS 
£3,467.2,313),/) 

INTX TINX TNX1 TNX2 TNX3', 

Pie TINY ONY] INY2 TNY3') 

10) NSTEP, INTX,TINX, INX, INTY, TINY, TINY 
(I5,4F7.2 


AAaANY 


9 


WRITE(6,971 
9712 FORMAT (I4,2( 

WRITE(6,970 
2700 oe 


ee oe 
14 


7 
2 
7 
# 
7 
9710 FORMAT Z 
END IF 


? 


Tit oelowine SEGINNING OF THE VGCMEAT CYCLE 


DO 4000 Je Poor 
NSTEP= 


AAAAN 


2en0 STORAGE FOR THIS TIMESTEP 


STNA=0 .0 
STNY=0 .0 
DO 1015 I=1,NTYPX 
ee} =0.0 
SBBB(I)=0.0 
1015 CONTINUE 
DO 1020 J=1,NTYPY 
SCCC(J)=0.0 
SDDD(J)=0.0 
1020 CONTINUE 


CALCULATE THE LOSSES DUE TO AREA FIRE FOR THIS TIMESTEP 


DO 1110 I=1,NTYPX 
DO 1105 J=1,NTYPY 
Sie) cee) {tanA(s, 1) ATNY(I)) 
sccc(s)=sccc(J)+(CcC(l , J )*TNX(I 
1105 CONTINUE 
1110 CONTINUE 
DO 1115 I=1,N 
AREX(1)= CINX(T) *SAAA(Z)) 
1115 CONTINUE 
DO #120 J=1, 
ARLY(J)= =TNY (3) *SCCC(3) 
21120 CONTINUE 


c CALCULATE THE LOSSES DUE TO AIMED FIRE PER TIMESTEP 


DO 1130 I=1,NTYPX 
DOVLIZ25 J=1 /NIYPY 
Sia jana habe Doo Beay Haney} 
SDDD(J)=SDDD(J)+(DDD(I,J)*TNX(I 
1125 CONTINUE 
1130 CONTINUE 
DOMES o I-! ,NIYPX 


AAMY 


AAD 


i3 


AMLX(1)=(INX(1)/TINX) *SBBB (I) 
1135 CONTINUE 
DO 1140 J=1 
pce (InY (3) /TINY) *SDDD(3) 
1140 CONTINUE 


c 
EC CALCULATE THE SELF LOSSES PER TIMESTEP 
e 


DO 1142 I=1,NTYPX 
SLX(I)=TNX(I) *UUU(I) 
1142 CONTINUE 
DON 1i4sea—i 
SLY(J) STNY (3) AVVV(S) 
1145 CONTINUE 


EC 
ce CALCULATE THE TOTAL KILLED BY ALL MECHANISMS 
C 


DO 1150 t=) ee 
TNDX(I)=-SLX(T) -ARLX (I) -AMLX(I)+NRRR(I) 
1150 Gentiiiia 
DO 1155 J—l Nive 
TNDY(J)=-SLY(3)- ARLY (J) -AMLY(J)+NSSS (J) 
1155 CONTINUE 


E 
: CALCULATE THE FORCE LEVELS AFTER THIS TIMESTEP 


DO 1160 I=1,NTYPX 
TNX(I)=TNX(1I)+TNDX(I) 
1160 CONTINUE 
DO 1165 J=1,NTYPY 
TNY(J)=TNY(J)+TNDY(J) 
1165 CONTINUE 


c 
C CALCUTE THE TOTAL X,Y FORCE LEVELS 
Cc 


DO 1170 I=1,NTYPX 
STNX=STNX+TNX (I) 
1170 CONTINUE 
TINX=STNX 
DO 1175 J=1,NTYPY 
STNY=STNY+TNY (J) 
1175 CONTINUE 


TINY=STNY 
C 
Cc 
C PRINT STATUS FOR EACH TIME STEP 
E 
WRITE(S 6710) NSTEP, INTX, TINX, TNX, INTY, TTNY, TNY 
6710 FORMAT(I4,2(15,4F7.2),/) 
c 
Cc 
C PRINT VALUES TO FILES 8,9,10,11 FOR PLOTTING 
Ec 
C 
G 


WRITE 8.8100} N, TINX , TINY 
WRITE(9,8100) N LINK (t) TNY (1) 
WRITE(10, 8100) N, TNX 2}, INE (2) 
WRITE(11,8100) N,TNX(3), TNY(3 

8100 FORMAT(13,2(F8.2 MQ. 
TERMINATE IF FORCES ARE ATTRITTED BELOW SPECIFIED LEVEL 
IF(TINX.LE.XSTOP .OR.TINY.LE.YSTOP) GO TO 5000 

RETURN AND CALULATE FORCE LEVELS FOR ANOTHER TIMESTEP 


4000 CONTINUE 


AAA AANM 
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AAADN 


5000 CONTINUE 


EXPLAIN CAUSE OF TERMINATION 


Ta Cte XSTOP) THEN 
E (6,9620) NSTEP,TINX, INTX 


ENDI r 
IF(TINY.LT.YSTOP) THEN 
WRITE (6,9630) NSTEP,TINY,INTY 


ENDIF 
IF (NSTEP .EQ. NSTOP) THEN 
WRITE (6, 9640) NSTOP 


ENDIF 
9620 FORMAT(' STOPPED IN STEP ',15,' BY FORCE lo), 
X ' CURRENT X STRENGTH =! amen INITIAL X SIRENGIH’ ,15,// ) 
E 


9630 FORMAT(' STOPPED IN STEP ',I5, ' BY FORC Ney, 
X ' CURRENT Y STRENGTH =' Bec). INITIAL WeoTRENGiH',15,//) 
9640 ae STOPPED BY PROGRAM AT NSTOP = BS) 


is 


APPENDIX E 
FORCE LEVEL TRAJECTORIES 


This Appendix is designed to provide the reader with a complete set of all of the 
results of this experiment. The input data sets for each model are presented in tables 


and then each of the force level trajectories is shown graphically. 


76 


a) 
(1) 
a 
O 
= 
yn 
oi 
ov 
O 
ee 
a om 
nD 
a 
= 
0 
a 
a) 
A 
0. 
(1) 
= 
QO. 
> 
< 
Sf) 


Ol NID Nix celocl loc2 loc3 U, Uy UZ Ty Ty T.3 
rel N ey? 1 loc! loc2 loc3 V1 V9 V3 8) S5 83 


TABLE 2 
SAMPLE INPUT DATA SET FOR M/W MODEL 
RaNIYeENIYS Vipv> v3 S] 55 53 


eseOr (2x2) 1 UNSTABLE ROOT § 
TEST OF (2X2), 1 UNSTABLE ROOT §$ 


TX1 NTX2 NTX3 vy Uy ug ry 1 13 


ke mene E MOIMeEHenen ONeN 
Fig D OW asco SAIN 
ee ake ato) Ooo *+ OOO gad - a © he ew ae ©) ay 

Mm ARAN AIAN AAA 
\or™ ey ey ecco cy 2 0,0 cacs 
~ A) ~-L — —— et oe ey — — 
= pt ep ON ONY TINO SION Sy NM aeennen FF MAM anon 
SL Amwe OOO 7.0000 a OAL. ODVVOOO LAOOOBSSH 


N 





di 


apie. S 
INPUT DATA SET FOR C/S MODEL CASE 1 


2 TEST OF (1X1), NO ROOTS 
1.0 1.0 1.0 0.040 0.000 0.000 28 0O 0 


().00000.0.00000 0.00000 
P00) OO Siro 

0.01000 0.00000 0.00000 
0.00000 0.00000 0.00000 
0.00000 0.00000 0.00000 
0.00020 0.00000 0.00000 
0.00000 0.00000 0.00000 
0.00000 0.00000 0.00000 

345215789.D0 45635761.D0 89342761.D0 


1.0 10-0010 0°000' 0:0CORaioS tan 


TABLE 4 
INPUT DATA SET FOR M/W MODEL CASE 1 


ou gs” TEST OF (1X1), NO ROOTS 


600 0 O 1.0 1.0 1.0 0.040 0.000 0.000 28 O Q 
Q. 96000 " QOQ000 OQ. DRGKG 


0 1.0 1.0 0.010 01000 0.000) 16 0 ae 
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TOTAL FORCE LEVEL TRAJECTORY 
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TOTAL X FORCE LEVEL 


CASE 1 (1X1) NO ROOTS 


Figure E.1 Total Force Level Trajectory For Case One. 
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X FORCE LEVEL TRAJECTORY 
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X FORCE LEVEL 


350 


° |....: Aan ae: LEGEND 
Sy} 2 i: | a C/8 MODEL 


200 





100 150 


012 3 4 8&8 6 7 8 9 10 11 12 13 14 15 16 17 16 19 20 
TIME 
CASE 1 (1X1) NO ROOTS 


Figure E.2. X Force Level Trajectory Over Time For Case One. 
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Y FORCE LEVEL TRAJECTORY 


700 


650 


LEGEND 
o C/S MODEL 


350 400 


Y FORCE LEVEL 
300 





250 


~~ 9 £2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18°19 20 
TIME 
CASE 1 (1X1) NO ROOTS 


Figure E.3 Y Force Level Trajectory Over Time For Case One. 
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TABLE 5 
INPUT DATA SET FOR C/S MODEL CASE 2 


De aiaik TEST OF (1X1), 1 STABLE ROOT AT 200,200$ 


1 600 O 0O 1.0 1.0 1.0 0.040 0.000 0.000 18 9O 0 
0.01000 0.00000 0.00000 ~ 
0.00000 0.00000 0.00000 
0.00000 0.00000 0.00000 

. .QQ000 O 


0.00000_0.00000 0.00000 
1 500 O O 1.0 1.0 1.0 0.060 0.000 0.000 22 0 QO 
0.01000 0.00000 0.00000 
0.00000 0.00000 0.00000 
0.00000 0.00000 0.00000 
0.00020 0.00000 0.00000 
0.00000 0.00000 0.00000 
0.00000 0.00000 0.00000 
345215789.D0  45635761.D0 89342761.D0 


TABLE 6 
INPUT DATA SET FOR M/W MODEL CASE 2 


01/05/87 TEST OF (1X1), | STABLE ROOT AT 200,200$ 


1 600. 0. Q. 0.040 0.000 0.000 18 O 9O 
00877 0.00000 0.00000 
00000 0.00000 0.00000 
0 0.00000 0. 


SS999 SSO9999 


0.00000 0.00000 0.00000 
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Figure E.4 Total Force Level Trajectory For Case Two. 
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TIME 
CASE 2 (1X1) 1 STABLE ROOT 


Figure E.5 X Force Level Trajectory Over Time For Case Two. 
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LEGEND 
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0 5 10 15 20 25 30 35 40 45 60 55 60 65 70 75 80 85 90 95100 
TIME 
CASE 2 (1X1) 1 STABLE ROOT 


Figure E.6 Y Force Level Trajectory Over Time For Case Two. 
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TABBEs] 
INPUT DATA SET FOR C/S MODEL CASE 3 


OL ge T TEST OF (1X1), 1 UNSTABLE ROOT AT 200,200$ 


1 500 0 0 1.0 10 10) ODIO" 000050 (0G am” 
0.06000 0.00000 0.00000 
0.00000 0.00000 0.00000 
0.00000 0.00000 none 


Q 1.0 1.0 0.010 0.000 0.000 18 O 90 


0.00020 0.00000 0.00000 
0.00000 0.00000 0.00000 
0.00000 0.00000 0.00000 
345215789.D0 45635761.D0 89342761.D0 


TABLE 8 
INPUT DATA SET FOR M/W MODEL CASE 3 


O1, eget TEST OF (1X1), 1 UNSTABLE ROOT AT 200,2008 


1 300. 0. Q. 0.010 0.000 0.000 22 O QO 
0.05 129 0.00000 9.00000 
0.00000 0.00000 0.00000 
Q. 00000 5 ORE K OQOQO00 
0.000 ae QO 
0.00000 Q 
0.00000 0, 00000 0. 00000 
1600. 0. Q. 0.010 0.000 0.000 18 QO 0 
0.03460 0.90000 0.00000 
0.00000 0.00000 0.00000 
0.00000 0.00000 0.00000 
0.00017 0.00000 0.00000 
0.00000 9.00000 0.00000 
0.00000 0.00000 0.00000 
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CASE 3 (1X1) 1 UNSTABLE ROOT 


Figure E.7 Total Force Level Trajectory For Case Three. 
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Figure E.8 X Force Level Trajectory Over Time For Case Three. 
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Y Force Level Trajectory Over Time For Case Three. 
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TABLE 9 
INPUT DATA SET FOR C/S MODEL CASE 4 


01/05/87 TEST OF (1X1), ROOTS AT 200,200(U) & 366,100(S)s 
0 1.0 1.0 1.0 0.040 0.000 0.000 28 0 0 


1 600 0 
0.06000 0.00000 0.00000 
0.00000 0.00000 0.00000 


000 0. 00000 0.00000 
1 600 0 10 1.0 1.0 OOLORUOO0. 0000 51a 080 
0.01000 a 00000 Q. 001000 


OO( 
345215789.D0 45635761.D0 89342761.D0 


TABLE 10 
INPUT DATA SET FOR M/W MODEL CASE 4 


QO. 0.040 0.000 0.000 28 9 0 


; 
.00000 0.00000 
. 0. 9. 0.010 0.000 0.000 12 0 90 


SS9999 SOS9S9 
Sea] COOOOO 
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Figure E.10 Total Force Level Trajectory For Case Four. 
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Figure E.11 X Force Level Trajectory Over Time For Case Four. 
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CASE 4 (1X1) 2 ROOTS 


Figure E.12 Y Force Level Trajectory Over Time For Case Four. 
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TABLE 11 
INPUT DATA SEI FOR. C/S3MIC DEL GASES 


Vary TEST OF (1X2), CASE 5 
000 0 1.0 1.0 1.0 0.030 0.000 0.000 38 O 0 


1.0 1.0 0.050 0.100 0.000 44 31 0 


SS9999 SOD999 


00 0.00000 
~ 345215789.D0 45635761.D0 89342761.D0 


TABLE 12 
INPUT DATA SET FOR M/W MODEL CASE 5 


02/02/87 TEST OF (1X2), CASE 5 
0, 0.030 0.000 0.000 38 0 0 


Q 


SS SSE RSS SSS 
SOOSO9SD OO9D90C— 
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Figure E.14 X Force Level Trajectory Over Time For Case Five. 
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Figure E.15 Y Force Level Trajectory Over Time For Case Five. 
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Figure E.16 X1 Force Level Trajectory Over Time For Case Five. 
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Figure E.17 Y1 Force Level Trajectory Over Time For Case Five. 


oo 


Y2 FORCE LEVEL 


Y2 FORCE LEVEL TRAJECTORY 


CCeceegeoes .- rre@reree gs Te eng se hee eg fete gy 0 Ff eget Seu ae ore ry ees wislce ne Cte cD (Meets 


950 


850 


ote Veccece® = e=Vege s blelr- cecil ele aS cielclelel\ieksinrers Nintsteiens 


eeeeetoee eesti ns C8evae 


es =H eee eeSeee cereogede ereege: 





650 


550 


F eeowere ee SMG OT CeCe SF Hw Hn Bee H Ct eC Cem uie ene S65 = e10\6 6:8 0.016 sieiel 6) 6)eisiciesieis) cevere'alisie eleiatslsieloiie Temeteretererenatene 


OQ ers }erete Por tog te eens eter eg cer g tee e tienes yee nee ps ct eg ews ceg On sr agecesssrenscgecscsges scp ences «cement 


S SMe a eeedeccn 5 Cee gsc eect Sectesesecn tis  selivisiseok el cieleieles Slelcvss e ¥cie i srebieteneretelaimeion SictelorieliiclcTainie Siena bias 
On eee SSS ESD OS 078 € 5 218 (oho ReReecnoner | me Sheil .7 e206) 0le1e 0) he) 0h mmmenens 


350 


eee eee eo teoeeeetee ese i s-e06 


150 





ee . ; : : 
0 & 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 
TIME 


0 


CASE 5 (1X2) UNSTABLE ROOT 


Figure E.18 Y2 Force Level Trajectory Over Time For Case Five. 
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TABLE 13 
INPUT DATA SET FOR C/S MODEL CASE 6 


Soe TEST OF (2X2), CASE 6, NO ROOTS 
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TABLE 14 
INPUT DATA SET FOR M/W MODEL CASE 6 
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Figure E.19 Total Force Level Trajectory For Case Six. 
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Figure E.20 X Force Level Trajectory Over Time For Case Six. 
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Figure E.21 Y Force Level Trajectory Over Time For Case Six. 
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Figure E.22 X1 Force Level Trajectory Over Time For Case Six. 
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Figure E.23 X2 Force Level Trajectory Over Time For Case Six. 
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Figure E.24 Y1 Force Level Trajectory Over Time For Case Six. 
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TABLE 15 
INPUT DATA SET FOR C/S MODEL CASE 7 


O1/31/87 TEST OF (2X2), CASE 7, 2 ROOTS 
: Pi meOomvOlO 0.020 0.000 11 16 0 
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TABLE 16 


INPUT DATA SET FOR M/W MODEL CASE 7 


01/31/87 TEST OF (2X2), CASE 7, 2 ROOTS 
y ee 0.020 0.000 11 16 Q 
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400. y 0.005 0.010 0.000 9 12 QO 
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Figure E.26 Total Force Level Trajectory For Case Seven. 
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Figure E.27 X Force Level Trajectory Over Time For Case Seven. 
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Figure E.28 Y Force Level Trajectory Over Time For Case Seven. 
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Figure E.29 XI Force Level Trajectory Over Time For Case Seven. 
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Figure E.30 X2 Force Level Trajectory Over Time For Case Seven. 
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Figure E.31 Yl Force Levél Trajectory Over Time For Case Seven. 
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Figure E.32 Y2 Force Level Trajectory Over Time For Case Seven. 
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TABLE 17 
INPUT DATA SET FOR C/S MODEL CASE 8 


PO! TEST OF (2X2), CASE 8; 1 STABLE ROOT 
emer O 100 0120 0.000 Sl 85 9G 


10 1.0 0.100 0.150 0.000 71 76 O 
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INPUT DATA SET FOR M/W MODEL CASE 8 
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Figure E.33 Total Force Level Trajectory For Case Eight. 


118 


700 


650 


X FORCE LEVEL 
450 500 550 


400 


350 


X FORCE LEVEL TRAJECTORY 


=e es @ spee wm pe en Ss oe Ow wrens 





t 
' 
a) 


iw 
— 
a 


es 


LEGEND 





0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95100 
TIME 
CASE 8 (2X2) 1 STABLE ROOT 


Figure E.34 X Force Level Trajectory Over Time For Case Eight. 
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Figure E.35  Y Force Level Trajectory Over Time For Case Eight. 
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Figure E.36 XI Force Level Trajectory Over Time For Case Eight. 
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Figure E.37 X2 Force Level Trajectory Over Time For Case Eight. 
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Figure E.38 YI Force Level Trajectory Over Time For Case Eight. 
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Figure E.39 Y2 Force Level Trajectory Over Time For Case Eight. 


124 


TABLE 19 
INPUT DATA SET FOR C/S MODEL CASE 9 


on TEST OF (2X2), CASE 9, 1 UNSTABLE 
2 4000 1.0 1.0 1.0 0.050 0.030 0.000 55 60 0 
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TABLE 20 
INPUT DATA SET FOR M/W MODEL CASE 9 
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Figure E.40 Total Force Level Trajectory For Case Nine. 
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X Force Level Trajectory Over Time For Case Nine. 
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Figure E.42 Y Force Level Trajectory Over Time For Case Nine. 
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Figure E.43 X1 Force Level Trajectory Over Time For Case Nine. 
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Figure E.44 X2 Force Level Trajectory Over Time For Case Nine. 
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Figure E.45 Yl Force Level Trajectory Over Time For Case Nine. 
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Figure E.46 Y2 Force Level Trajectory Over Time For Case Nine. 
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TABLE 21 
INPUT DATA SET FOR C/S MODEL CASE 10 


fee OU PESTOr (3x5) 


Pee007150 175 10 1-0 1.0 0.005 0.010 0.020 4 4 II 
0.06000 0.02000 0.02000 
0.00000 0.00000 0.00000 
0.03000 0.01000 0.02000 
0.00000 0.00000 0.00000 
0.00000 0.00030 0.00020 


5 1.0 1.0 1.0 0.003 0.020 0.030 4 5 14 
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TABLE 22 
INPUT DATA SET FOR M/W MODEL CASE 10 
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Figure E.47 Total Force Level Trajectory For Case Ten. 
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Figure E.48 X Force Level Trajectory Over Time For Case Ten. 
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Figure E.49 Y Force Level Trajectory Over Time For Case Ten. 
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Figure E.50 X1 Force Level Trajectory Over Time For Case Ten. 
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Figure E.51 X2 Force Level Trajectory Over Time For Case Ten. 
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Figure E.53 Y1 Force Level Trajectory Over Time For Case Ten. 
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Figure E.54 Y2 Force Level Trajectory Over Time For Case Ten. 
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Figure E.55 Y3 Force Level Trajectory Over Time For Case Ten. 
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